THE BULLETIN OF 
Mathematical 


BIOPHYSICS 


SEPTEMBER 1956 


What Type of Empirically Verifiable Predictions Can Topological 


Biology Make? —N. Rashevsky- - ------+--+-+-+-++-+-.-- 173 
A Neural Mechanism for Adjustment to Optimal Conditions, With Pos- 
sible Reference to Visual Accommodation—N. Rashevsky- - - - - 189 


Diffusion and Simultaneous Chemical Reactions: IIl. The Degree of 
Localization Achieved in Cytochemical Staining Procedures—D. _ 
G. O'Sullivan - - = - - - +--+ - ee ee eee ee ee ee 199 


_A Probabilistic Approach to Some Problems in Blood-Tissue Exchange 


sRobert Macey® = \- = 050 = eo) es ww eee tw na aie 205 ees ; 


Parasitism and Symbiosis in an N-Person Non-Constant-Sum Contin- 29 
_ vous Game—Caxton Foster and Anatol Rapoport - - - - ---- - 219 


The Geometrization of Biology: A Correction—N. Rashevsky - - - - 233 
On The Information Content of Gravha: Compound cel Different 


States for Each Point—Ernesto Trucco - - - -----+--+--- 237 
On The Diffusion ot lons in Membranes—Joseph D. Harris- - - - - - 255 
Errata aes tira eatin ewe eo fa latieiere bite a lie bose) mica fe 263 


THE UNIVERSITY OF CHICAGO PRESS - CHICAGO 
VOLUME 18 arts 7 Bees Siege ; NUMBER 3 


The Bulletin of 
MATHEMATICAL BIOPHYSICS 


Eprror: 


N. RASHEVSKY, UNIVERSITY OF CHICAGO 
AssociaTE EpiTors: 
H. D. LANDAHL, UNIVERSITY OF CHICAGO 
ANATOL RAPOPORT, UNIVERSITY OF MICHIGAN 


The BULLETIN is devoted to publications of research in Mathematical 
Biology, as described on the inside back cover. 


Tue BuLLETIN is published by the University of Chicago at the University of 
Chicago Press, 5750 Ellis Avenue, Chicago 37, Illinois, quarterly, in March, June, 
September, December. {The subscription price is $10.00 per volume, the price of 
single copies is $3.00. Orders for service of less than a full volume will be charged at 
the single-copy rate. Postage is charged extra as follows: for Canada, 20 cents per 
volume (total, $10.20); for all other countries in the Postal Union, 50 cents per vol- 
ume (total, $10.50). YPatrons are requested to make all remittances payable to the 
University of Chicago Press in postal or express money orders or bank drafts. 


THE FOLLOWING is an authorized agent: 


For the British Commonwealth, except North America, aig and Australasia: 
The Cambridge University Press, Bentley House, 200 Euston Road, London, N.W. 1. 
Prices of yearly subscriptions and of single copies may be had on application. 


CLAIMS FOR MISSING NUMBERS should be made within the month following the 
regular month of publication. The publishers expect to supply missing numbers free 
only when losses have been sustained in transit, and when the reserve stock will 
permit. a . 

BUSINESS CORRESPONDENCE should be addressed to the University of ee 
Press, Chicago 37, Ill. 


COMMUNICATIONS FOR THE EDITOR and manuscripts should be addressed to N. 
_ Rashevsky, Editorial Office of The Bulletin ‘e Mathematical Biophysics, 5741 Drexel 
Avenue, Chicago 37, Il. “ 


NOTICE TO SUBSCRIBERS 
If you change your address, please notify us and ina local postmaster immediately. 
[Copyright 1956 by the University of Chicago] © 


Permission to reproduce for purely scientific or scholarly purposes ri material | 
published in this journal will be given upon request. 


Re-entered as second-class matter April 6, 1956, at the Post Office at 
Lancaster, Pennsylvania, under the Act of March 3, 1879, 


BULLETIN OF 
MATHEMATIC AT. BIOPHYSICS 
VOLUME 18, 1956 


WHAT TYPE OF EMPIRICALLY VERIFIABLE PREDICTIONS 
CAN TOPOLOGICAL BIOLOGY MAKE? 


N. RASHEVSKY 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


It has been shown previously (Bull. Math. Biophysics, 17, 207-218, 
1955) that within the framework of the general principle of histepalogical 
mapping a number of different specific hypotheses can be made about the 
topological structure of the primordial organism and the transformation 
which lead to the description of higher organisms. Each set of such hy- 
potheses leads to different conclusions and different predictions of spe- 
cific biological relations. In the present paper it is shown that while the 
specific predictions depend on the specific assumptions within the frame- 
work of the general principle of biotopological mapping, the type of the 
predicted relation is characteristic of the general topological approach. 
Relations of this kind would hardly ever be considered in the usual metric 
approach to mathematical biology. As specific examples it is shown how 

_topological considerations lead to the problem of the total number of pos- 
sible organisms, the problem of relations between number of organs and 
number of cell types in an organism, the problem of the number of cell 
types in different organs, a possible approach to the problem of the num- 
ber of different endocrine secretions in different organisms, and other 
similar problems. All this shows the predictive value of topological 


biology. 


In previous papers (Rashevsky, 1954; 1955a,b,c; 1956a,b; here- 
jnafter referred to as IJ, II, Il, IV, V, and VI) we outlined a topo- 
logical approach to biology based on the principle of biotopological 
mapping (I, pp. 324-25) which is inferred from some generally known 
facts of biology. It has been pointed out in I that the metric ap- 
proach, which has been used hitherto in mathematical biology and 
which has proven its predictive value, has serious limitations. It 
leaves entirely out of consideration the numerous non-metric or re- 
lational aspects of biology, which are in many cases just as im- 
portant as the metric ones. Topology, as the mathematics of rela- 
tions rather than of quantities, is much more adapted to the de- 
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scription of the relational aspects of biology. This better adapta- 
tion, however, is not sufficient. It must also be shown that topo- 
logical biology can lead to empirically verifiable predictions of 
such a type as cannot or are not likely to be made by the metric ap- 
proach. The final proof must be made by actually exhibiting some 
predictions which are verified empirically. This goal may still be 
far removed, awaiting a much more complete development of topo- 
logical biology. What can be done at this early stage, however, is 
to show what type of predictions the different systems of abstract 
biologies (IV) can make. It is less important now that those pre- 
dictions should actually agree with empirical data than that they 
should have the following two properties: (a) They should indicate 
new ways for experiment and observation, and (b), if a prediction 
turns out to disagree with empirical data, this disagreement must 
indicate the way in which some basic assumptions of the theory 
must be changed in order to obtain agreement. To give a few ex- 
amples of such types of possible predictions is the purpose of this 
note. 

Four verifiable conclusions have already been drawn even at 
this early stage of development: one is the total number of different 
possible genera of organisms that can exist (I, pp. 341-42); the 
other is that the more differentiated an organism, the less its power 
of regenerating lost organs and functions (Il, pp. 124-25); the 
third, that more complex organisms are more adaptable to losses of 
some parts or functions, as well as to changes in the environments 
(IV, pp. 211 and 213); the fourth is a relation between the total 
number of different cell types and the number of organs of an or- 
ganism [IV, equation (1)]. 

Of the four conclusions, the second and third represent two 
known general biological laws and thus are empirically verified. It 
has been pointed out in IV that the agreement is likely to be to 
some extent accidental and that at this stage too much value should 
not be put on the actual agreement. We have here two laws, how- 
ever, which are semi-quantitative in nature and which, for reasons 
discussed in I, are not likely to be derived from the metric approach 
to biology. 

As has been shown in IV, a choice of a given primordial graph 
and of a given transformation determines a system of abstract 
biology. Of all possible systems of abstract biology one comes 
nearest to the description of actual biological phenomena, just as 
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amongst the different systems of abstract geometries one is the 
most useful in describing our physical universe. As long as there 
is no general theory of systems of abstract biology, the problem of 
choosing the proper one by a sort of mere guess is practically in- 
superable. However, if we develop various systems of abstract 
biology, deriving for each system a set of theorems, then by com- 
paring the sets of theorems thus derived with the set of empirical 
data, we shall choose that system as the best which shows the 
best agreement between theory and observation. Thus at first a 
purely abstract study of different conceivable systems of abstract 
biology is unavoidable. 

We have seen in IV that some properties of a system of abstract 
biology are determined only by the choice of the transformation, re- 
gardless of the choice of the primordial graph. Others are de- 
termined by both. Systems which depend only on the choice of the 
primordial graph are of little interest since they deal only with the 
properties of the primordial organism. In VI we have studied alto- 
gether 14 different transformations T™ (u=1,..., 14). Two dif- 
ferent primordial graphs, one represented in Figure 1 of I, the 
other on Figure 4 of VI, have been hitherto considered. We shall 
denote them respectively by P, and P,. We shall denote a sys- 
tem determined by a given transformation T and a given pri- 
mordial P, by (TOP. y, indicating, if necessary, in additional 
parenthesis the range of variation of the indices 7 and &. Thus we 
may have r@R,) (i= 1,...,8; K=1, 2). The class of systems 
which is determined by a given T® and any P, shall be denoted by 
(TX). Thus if a theorem holds for (TX), it means that it holds 
for all systems based on T® regardless of the choice of P,. Thus 
theorem 2 of III (or, as we shall denote it, theorem III.2) with the 
correction made by E. Trucco (1956) holds for (TX); so do theo- 
rems IV.1 and IV.2. 

I. Let us now dwell on the other two conclusions. As remarked 
in I, the number of organisms, or rather genera, estimated theoreti- 
cally for (aes) is about 1000 times greater than the actual num- 
ber observed. While it is quite possible that only a very small 
fraction of all the possibilities have hitherto been realized in na- 
ture, it seems also very likely that the topological transformation, 
which we denoted by T in I and by T in VI, is too general for our 
purposes. Moreover the number of organisms obtained from 7“ can 
only be estimated very roughly as no closed expression for that 
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number can be derived. The total number is determined basically 
by two factors. First by the number of different ways in which n 
specializable biological functions (I) can be distributed between m 
different cell classes or, as we may call them somewhat loosely, 
tissues. This number FR” is given (II) by 


so ie eee ge pee 
a joie ‘Gap pe s 


Since for a fixed n, m can vary from 1 to n, therefore, the total 
number of possibilities is given by 


n 
S,= ). RR. (2) 
m=l 

In topological terms; R” represents the number of ways in which 
nm specializable points can be distributed among the m component 
primordial graphs, while S, represents the total number of possi- 
bilities for a given n.- 

Within each component primordial graph we do have a finite num- 
ber of possible distributions described by Te (I, p. 334). The 
number gw of such distributions is given (II, p. 123) by 

t 


Cal n—-1 
Q. P = n,=~1 ’ (3) 
u t 
the notations being the same as before. 
There is, however, no way of dividing either S, or RF” into 
‘classes characterized by the same gk In fact in each component 


primordial graph, for any m there is in general a different number 
Q"~', Hence we cannot merely multiply S$, by Qn this would 


make no sense because of the variable n,. The only way of esti- 
mating the total number of possibilities is to take the value of s, 
and to multiply it by some very roughly estimated average of ape : 
The value of S, increases very rapidly with n as seen from the fol- 
lowing: S, = 5; 8S, = 15; S, = 52; S, = 203; S, = 877. 

The estimates mentioned on page 341 of I were made on the basis 
of the primordial graph shown in Figure 1, page 326 of I, for which 
m=17 (I, p. 341). Considering that the actual differentiation of 
animal and vegetative functions occurs always in distinct groups, 


RT Tbe ee 
. x 
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we first considered S19 Which corresponds to the 12 vegetative 
functions shown in Figure 1 of I, and then S;, corresponding to the 
7 animal functions, eS My C, 1, M,, S,, and M,- For each of the 
S,) possibilities of distributing 10 points amongst different tissues 
(m=1,..., 10), there are S, possibilities of distributions of the 7 
animal points. Altogether S, x S,, ~ 10°. 


As for an average value of Gee we may argue like this: The 

i 
value of n, varies from 1 to n. Since’we subdivided the points in 
two groups, one with n = 7, the other with n = 10, we may consider 


as a plausible average for QO? the value Os. = : = 126 and 
t 


6 
for ge the value ase = (*) = 15; a total average would be of the 
i 


order of fifty. Then we find for the total number of transformed 
graphs 


50 x-10° =.5 10°. 


The figure of 10° arrived at in I, page 341, is obtained by con- 
sidering the following biological argument. Some related biological 
functions do not differentiate independently. There would be no 
physical advantage for an animal to possess a highly developed 
eye which can perceive small details if its locomotion mechanism 
is not developed sufficiently enough to respond to those details. 
Similarly, the development of other parts of the nervous system is 
correlated with the development of the motor apparatus. Thus 
specialization occurs actually in ‘‘blocks’’ of points. 

This phenomenon has the same effect as a reduction of n. Con- 
sidering in Figure 1 of I that the points S, and M,, as well as the 
points F’., M C, and M,, develop in a correlated manner, the num- 
ber of possibilities for the animal function is reduced from S, = 877 
to S, = 5, with a corresponding reduction of the final result to 10°. 
The whole argument is neither accurate nor at all convincing. 

It is readily seen that the order of magnitude thus estimated is 
the same for all the 14 transformations T™ (u=1,2,..., 14) so 
that all the (T™P,) (w=1,..-, 14; ¢= 1,2) lead to an extremely 
high value of N, the total possible number of biotopologically 
distinct organisms. 

As remarked above we cannot strictly affirm that such a high 
value of N contradicts the empirical data, Yet it is highly probable 
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that it is so. Now suppose, just as an illustration, that in some 
way it has been established that the actual number N is much less 
than the one estimated above. Now, should we modify our system 
of abstract biology to bring the theoretical value of N into better 
agreement with the observed one? 

In all systems (7™X)(u=1,..., 14) the number n of speciali- 
zable points is fixed. Moreover, for any m all n points are special- 
jzed in one of the component primordials, according to Foes (I, VI). 
We may now restrict the transformations Te) (w=1,...; 14) by as- 
suming that of the n specializable points only m < n are actually 
specialized for a given m. This means that for any given choice of 
m as parameter of the transformation some m of n specializable 
points are actually specialized, that is, each component primordial 
loses m—1 points, specializing only in one rather than in n, bio- 
logical functions. Theremaining n —m specializable points remain 
non-specialized, that is, they are not lost by any of the m com- 
ponent primordials. Each component primordial receives now only 
one specialized point, losing m-—1 others. The step Hod (u = 
1,..., 14) is now eliminated while the step Le becomes much 
simpler. Thus the one specialized point receives now m— 1 sub- 
sidiary points. 

The maximum amount of specialization is obtained when m =n, 
when all n specializable become actually specialized. The resid- 
ual graphs, attached to the subsidiary points should not contain 
any specializable points. To each transformation 7™ (aa Tees 
14) we now obtain a corresponding simplified transformation T™’. 
The number of possibilities for any T™” is easily calculated. For 


a given m there are a possibilities, namely, the number of pos- 


sible choices of m points out of n. Since m varies from 2 to n, 


therefore in (T™’X) (w=1,..., 14) the number AN of all possible 
biotopologically different organisms is 


Woes (7) =2"=n=1. (4) 


In P, we have n= 17 and therefore in (TPs (wettest 
we have N = 10°; in P, we have n= 22 and therefore in (TSP 
we have N = 10". 
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If we again consider the circumstance that some of the n points 
do not specialize independently, this again reduced n to what may 
be called effective n,. The value of nm, is to be substituted for n 
in equation (4). In P, this reduces nto n_ = 10, giving N ~ 10° for 
(T”P,) and N = 10° for (T™”P,), . 

What is the biological meaning of substituting T™* for T? In 
(TX) several biological functions specialize in one group, giving 
what we may call a polyfunctional organ, that is, an organ which 
performs several specialized biological functions. In (POO each 
organ performs only one specialized biological function. Inasmuch 
as we have such apparently polyfunctional organs as the skin, 
which performs the biological functions of sense of touch, sense of 
heat, and the biological function of elimination, it may seem that 
T™’- should be discarded at once in spite of some apparent ad- 
vantages. However, we must consider the following. In (T™X) all 
the specialized functions of a polyfunctional organ are lost by all 
other parts of the body. In (7™’X) an organ may still be poly- 
functional, but of all the biological functions performed by it only 
one is completely specialized and lost by all other parts of the 
body. The remaining ones may be performed also by other organs 
or parts of the body. Considering the fact that elimination is not 
a monopoly of the skin and that sensitivity to touch is preserved by 
some internal parts of the body, we see that not only can T™” not 
be summarily discarded but that, on the contrary, it may deserve 
more consideration. 

In any case the problem of the number of possible organisms, 
which appears in a natural way at the very beginning of topological 
biology, is not one that is likely to be treated successfully by the 
metric approach to mathematical biology. 

Il. We now turn our attention to equation (1) of IV which gives 
us a relation between the number 7 of different tissues, the number 
m of organs, and the number n of specializable biological func- 
tions. Equation IV.1 holds for (TX) (w=1,..., 7). In the pre- | 
ceding papers we used the word tissue to designate a connected 
group of cells of identical type. Though this definition of the word 
is pretty much in line with the usage in biological and medical 
sciences (Stedman, 1946, p. 1154) yet some histologists use the 
word in a much more restricted sense, speaking only of four types 
of tissues: nervous, muscular, epithelial, and connective. To avoid 
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misunderstanding we shall refer to the quality 7 in equation IV.1, 
which we repeat here, 


r=m+n(m-—1), (5) 


as number of cell types. 

Equation (5) above is verifiable empirically. It may even seem 
to be easily verifiable as 7, and m could be directly counted. For 
a given organism, we then could determine n. With this value we 
should obtain for all other organisms which are derived from the 
same primordial a relation between r and m, free of arbitrary param- 
eters. Since the number of primordials is either one or is very 
small the verification of (5) would seem to present no particular 
difficulties. 

Actually the situation is not quite so simple, because neither r 
nor m has ever been systematically counted for different organisms. 

Therefore there are no data now available in the literature which 
would enable us to verify (5) over a wide range of organisms. The 
method outlined above is, however, valid, and whether equation (5) 
proves to be correct or not, it definitely suggests a new avenue of 
experimental research. It is again a problem that seems quite for- 
eign to the metric approach. 

In (T™’°X) (w= 1,..., 7) equation (5) is simplified by making 
n=m, to 


r=m°, (6) 


that is, the number of cell types of an organism is equal to the 
square of the number of organs. Since neither (TX) nor (T™’X) 
are likely to correspond to biological reality, we should not expect 
either (5) or (6) to hold actually. 

From equation VI.19, it follows that in (T™X)(u = 8,..., 14) we 
have instead of (5) 


r=m+(n-n,)(m-1), (7) 


where n, is defined in VI. In the corresponding (T™’X)(u = 8 
14), we have instead of (6) 


geeeg 


2 
T= M -NpM+ Np - (8) 


A systematic observational search for a relation between r and m 
would be well worth while and would throw light on the biotopologi- 
cal relations in the organic world. 
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Il. An important class of specializable biological functions is 
formed by internal secretions which play a coordinating role in the 
over-all functioning of the organism. More and more hormones are 
being discovered, and a few dozens are known now in humans. 
How many hormones are there altogether in a given organism? This 
is a question which hardly has a place in metric mathematical 
biology, but to which topological biology may perhaps give an 
answer some day. To suggest a physico-chemical scheme for the 
chemo-regulatory processes which regulate all the hundreds of 
separate physiological functions of a complex organism is a hope- 
less task. To suggest a similar scheme for a relatively simple 
primordial organism may be quite feasible. We may, for example, 
assume that each of the specializable biological functions of a 
primordial organism is regulated by a special hormone or that one 
hormone coordinates the biological functions in a ‘‘block’’ of such 
functions which specialize in a correlated manner. Relatively 
simple interactions between the hormones of the primordial may be 
postulated. Once this is done, a proper transformation T™ will 
give us the number of subsidiary biological functions to a given 
specializable internal secretion, and those subsidiaries are to be 
interpreted also as additional internal secretions. Moreover T™ 
also gives us all the connections between the points which repre- 
sent the different additional internal secretions. As we have seen 
in I, a very simple partial of the primordial is, in general, trans- 
formed into a very complicated partial of the resulting graph. A 
guess at a correct partial graph for the internal secretions of a 
complex organism is a hopeless task. But such a correct complex 
graph may well be derived from a much simpler primordial which is 
much easier to postulate. 

Adding even only a few internal secretions to the primordial will 
increase its complexity and therefore increase the number N of 
possible organisms. With specialization in ‘‘blocks’’ we obtain, 
however, for (TEP) Gi = 1,..., 7) a value N ~ 10° (p. 179) which 
' is much too small. Hence there may well be room for the addition 
of more points to the primordial. 

If, out of n specializable points of the primordial, n, stand for 
internal secretions, then in a completely differentiated graph, in 
which m =n, the total number of points will be n” since each of 
the n points receives n subsidiaries, whereas the number of points 
which represent the internal secretions in the completely differenti- 


189 N. RASHEVSKY 


ated graph is equal to nn,. Hence the ratio of the number of in- 
ternal secretions to the total number of specialized organs is the 
same in the completely specialized graph as in the primordial. In- 
asmuch as we remarked above that even in humans there are poly- 
functional organs, whereas a completely differentiated organism 
could have only monofunctional organs, we see that no completely 
differentiated organism has yet been evolved. But in intermediate 
steps the ratio of the number of internal secretions to the total 
number of organs is not equal to n,/n. In fact, for the same value 
of m it may vary depending on the distribution of the specializable 
points between the m component primordials (Step POS I, p. 334). 
A simple expression for all the possibilities can be derived in 
(FOX) Set) 

In T™’ only m<n of the specializable points actually spe- 
cialize. The total number of possible choices of those m points 


out of nis @ , when all points specialize independently. Amongst 


the n points there are n, which represent the internal secretions. 
In how many ways, K, can we choose m points out of n in such 
a way that & < m of those points will be from the group of 
the n, points, which represent the internal secretions? There are 


fe) ways of choosing & points out of n,. For each choice there 


remains n—mn, points, not belonging to n,, from which we can 
choose m — & points to be added to & making a total of m. But there 


are fe eo ":) ways of making that latter choice. Hence the answer 


m—k 
Reet area): @ 


to our question is 

For each possibility in which & of the n, points are specialized, 

each of the & points gives rise to m — 1 subsidiaries or to a total 
of m points each. The total number of points generated by the & 
points is mk. Hence the total number of possibilities of generating 
mk points is given by (9). 

We now can compute the average number n, Of points which cor- 
respond to internal secretions for all organisms with a given m, as- 
suming equal probability of all combinations. To this end we mul- 
tiply (9) by mk and sum for all values of &, that is, from & = 1 to 
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kK=m,ifm<n,; and fromk=1tok=n, if m> n, Then we divide 
the result by the total number (") of possible combinations of n 


point by m. We thus find 


& Oo8) 


Ny 5 , form<n,, (10) 
(n) 
and 
y" G - z) (7) at 
< & m—k k 
Ny, , form>n,. (11) 


The value of 7, is a function of m, n, and n,. Thus (10) and 
(11) give us the average number of internal secretions in the class 
of organisms characterized by a given number m of different or- 
gans, for a system ease (w=1,..., 7) in which P is arbitrary 
except for containing n specializable points, of which n, < n repre- 
sent hormonal secretions. Here again we have a prediction of a 
type of relation which has never been looked for in biology. 

Corresponding relations for Cer) (w=1,..., 14) are much 
more complicated. 


IV. In (T“’X) (u=1,..., 7) there is only one specialized point 
per component primordial, and it has m~-—1 subsidiary points to 
each of which is attached a residual graph. Hence all organs are 
monofunctional and have the same number of residual graphs for 
each component primordial. In biological terms we have the 
following: 

Theorem 1. In (T™’X) (w=1,..., 7) the number of cell types in 
all organs is the same and is a function of the number of organs 
only. 

It does not seem likely that such a relation actually holds in 
nature, but it is certainly an empirically verifiable relation. 

The situation is different for (TX) (u=1,..., 7). Here if a 
component primordial contains n, specialized points, it contains a 
total of n-—n, subsidiary residual graphs, or, together with the 
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given component primordial, a total of n—1n,+ 1 partial graphs. 
Hence all component primordials, which contain the same number 
n, of specialized points, though the actual points may differ from 
one component primordial to another, give rise to the same number 
of partial primordial graphs. In biological language we have 
Theorem 2. In(T™X) (u=1,..., 7) all polyfunctional organs of 
an organism, which perform the same number of different biological 
functions, consist of the same number of different cell types. 


V. In IV, theorem 2, we have seen that in (TX) if two non- 
adjacent specializable points of the primordial are connected by a 
way (Ill, p. 114), then the same two specialized points are con- 
nected in the transformed graph by at least m ways. One biological 
consequence drawn from this theorem was that the more complex 
organisms are more adaptable to possible losses of some biological 
functions. There is, however, another important consequence of 
that theorem. Let m ways lead from a specialized point f, to a 
specialized point 7,. If the biological function represented by f, 
is in some manner impaired, this impairment will affect all biologi- 
cal functions to representative points of which there is a way from 
the point f;. But that means that more biological functions will be 
affected in a complex organism than in a simpler one, due to im- 
pairment of the same specialized biological function f;. Any af- 
fected biological function manifests itself either as an objective 
sign (Stedman, loc. cit., p. 1029) or, in humans, as a subjective 
symptom. Thus many signs or-symptoms will be produced by the 
impairment of a biological function. A complex of symptoms con- 
stitutes a syndrome (Stedman, loc. cit., p. 1105). We shall use here 
the word syndrome somewhat loosely as referring both to signs and 
symptoms. Therefore, from Theorem IV.2 there follows 

Theorem 3. In (TX) the complexity of a syndrome due to the 
impairment of a specialized biological function increases with in- 
creasing complexity of the organism. 


It is readily seen that theorem IV.2 and therefore also theorem 3 
hold in (TX)(u= 1,...,14) as well as in (T’X)(u=1,..., 14). 
This follows from considerations which lead to equations (21) and 
(22) of VI. 


We now shall prove a theorem which at the present stage of de- 
velopment is probably biologically irrelevant, but which illustrates 
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another interesting possible application of topological con- 
siderations. 

Let us consider the so-called fundamental sets (Grundmenge of 
Konig, 1936). A fundamental set is defined in the following way. 
Denote by II, the set of points which can be reached by a directed 
way (Rashevsky, 1955b, p. 114) from the point a. If II, is not a 
proper subset of any other Il,, where 6 is not included in IT, then 
IT, is called a fundamental set G,, with the point a as its source. 
If directed circles are present, then a fundamental set may have 
more than one source. In general, different fundamental sets of a 
graph have a non-empty intersection. Some of them may, however, 


FIGURE 1. 


be disjoined. For example in the graph of Figure 1,G, consists of 
the points 1, 2, 3, 4, 5, 6, 7; G, consists of 8, 3, 6, 7, 9, 10; 
G,,—of 11, 9, 10. G, and G,, are disjoined. G, and G, have a 
non-empty intersection consisting of points 3, 6, 7; G, and G,, 
have a non-empty intersection consisting of points 9, 10. 

If the sets G, and G, are disjoined, then no point of G, can be 
reached by a way from any point of G, (7 #7), and vice versa. 

Our theorem is stated as follows: 

The number of mutually disjoined fundamental sets, each having 
a specialized point as its only source, is an invariant of the 
transformation (pede 

Proof. Let the primordial graph contain & mutually disjoined 
fundamental sets G, (¢=1, 2,..., k) with corresponding sources 
S, (¢=1,2,...,%). The transformed graph T consists of m com- 
ponent primordials, which all have the points S, in common, and’ 
of a number n(m—1) of subsidiary residual graphs (Rashevsky, 
1954, p. 335). 
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In the primordial graph P let the fundamental set G, consist of 


7, points a, » a; orale a, s Denote the points which correspond to 
i 


t 
a, in any of the m component primordial graphs by a; (9 (g=1,2,. 
m). Then in 7” there are mr, points ae, a9 a, oe. ; ae Ge 1. 
r 
i 


Q,..., m), which can be reached by a directed way from S,. In 
addition some points in the subsidiary residual graphs which be- 
long to S, according to T,, and T, (Rashevsky, 1954, p. 334) may 
be reached from S,. 

a) We shall now prove that all those points together with the 
point S, itself form a fundamental set Be of T‘, Suppose this 
were not the case. Then there would exist at least one point A of 
T not belonging to Ge such that Ge would be a proper subset of 
Sie But since ait the points aed all the directed ways of T“ 
can be mapped on corresponding points and ways of the primordial 
P, therefore, the primordial would contain a point 4,, on which A 
is mapped, such that A, would not belong to G, and such that G, 
would be a proper subset of Ty . But then G, would not be a 


fundamental set, which is contrary to the assumption. 

Thus the mapping of a set ee , which is not a fundamental set 
in T™, onto P gives also a non-fundamental set in P. By a similar 
argument we prove that the mapping of any fundamental set of 7“? 
onto P gives a fundamental set in P. 

b) Next we shall prove that all fundamental sets GP are mutu- 
ally disjoined. If GP and GR (¢ #.l) were not digjeined then there 
would be at least one way in T“ which would lead from a point 

ae of Ge to some point ae of Ge But this way would map on 
a way from a, of G, to a, a G,. On the primordial P this would 
mean that G, and G, are ‘not disjoined, which is contrary to the 
assumption. 

From a and 6 it follows that to each G, (¢=1, 2,..., k) of P 
there corresponds a ee of T with the ae source S,. Hence 
there are at least k fundamental sets rar in T with Ee 
ing sources 8S @erty 2s. ce, &). 

c) Now we ae complete our proof by showing that there are not 
more than & disjoined fundamental sets Ge Such that each has a 
specialized point as its only source. Suppose there were a(k + 1)st 
fundamental set Ca with a specialized point 8,4, 8S source, and 
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that ce were disjoined with all other fundamental sets Gee . All 
the points and ways of oy will map on corresponding points and 
ways of P. Hence, contrary to the assumption, P would contain a 
(& + 1)st fundamental set, disjoined from all other G's. 

Thus our theorem is demonstrated. 

Let us now translate it into biological language. As remarked 
above, if the biological function represented by the point 8, is 
disturbed, this will result in a disturbance of all biological func- 
tions which are represented by the other points of G,. All these 
disturbances form a syndrome. We shall call two syndromes dis- 
joined if they do not have any symptoms in common. Our theorem 
then means biologically that the number of disjoined syndromes, 
due to disturbances of some specialized biological functions, is 
the same for all organisms. Inasmuch as the theorem has been 
proved as valid only for the transformation 7“, which is a biologi- 
cally implausible one, therefore no direct applications of the above 
are to be expected. Moreover, practically no one has ever bothered 
to study syndromes of any kind in lower animals. However, when 
it comes to the all important problem of extrapolating results of 
medical experimentation on higher animals to human beings, the 
study of relations of such type as the above may prove to be of 
great practical value. 


This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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A general mechanism is studied which adapts the value of some bio- 
logical parameter so as to maximize or to minimize the value of another 
parameter. In the process of visual accommodation, the thickness of the 
lens is determined by the maximum sharpness of the image on the retina. 
A neural mechanism which provides for such an optimum adjustment is de- 
scribed. A similar mechanism may play a role in such phenomena as 
hedonistic behavior where activity is adjusted to a maximum of satisfac- 
tion. The mechanism discussed here does not provide for an absolutely 
sharp retinal image, but leads to a conclusion that the sharpness of that 
image fluctuates within narrow limits which are determined by the param- 
eters of the systems. 


In the study of hedonistic behavior (Rashevsky, 1951, chap. iv) 
we suggested a neural mechanism which adapts the intensity of 
some activity of an organism to such a value which corresponds to 
the maximum value of satisfaction, when the latter is considered 
as a function of the intensity of the particular activity. Such situ- 
ations occur in many other fields of biology. A possible case, for 
example, is visual accommodation. The thickness of the lens, 
which determines its radii of curvature, is maintained by the ciliary 
muscle at such a value at which the image on the retina is the 
sharpest. The greatest sharpness may be said to be equivalent to 
the maximum contrast, C, or to the minimum size, d, of the retinal 
image, or to the maximum information which the retina can receive 
from the object. The information-theoretical point of view may in 
the long run turn out to be the most interesting one. We shall, 
however, choose here the contrast C. 

The exact nature of the quantity the maximum of which deter- 
mines accommodation is immaterial for the following considera- 
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tions. We shall denote it by C and not discuss it in any more de- 
tails. All that matters is that the thickness of the lens is adjusted 
so as to maximize C’. 

This quantity C is a function, C(r), of the thickness r of the 
lens and has one maximum. The thickness r remains stationary 
when dC/dr = 0; if dC/dr > 0, then r increases so that dr/dt > 0; 
when dC/dr < 0, then dr/dt <0. From a formal point of view, the 
problem is analogous to that of maximizing the satisfaction func- 
tion (Rashevsky, 1951, chap. iv). Therefore, formally we may set 


(loc. cit., p. 42) 
d dC 
ae (=), (1) 


where F(z) is a monotonically increasing function of 2, such that 
F(0) =.0. 

Again in a general theory of adjustment to optimal conditions 
the nature of 7 is irrelevent. If C denotes, for example, satisfac- 
tion and r the intensity of the activity which produces (C, then 
equation (1) still formally stands, although it is now interpreted as 
IN; l0C;<c2t (1951). Thus, mutatis mutandis, our discussion below 
may be applied to any biological phenomenon in which a variable r 
is adjusted to such a value as to maximize another variable C, 
considered as a function of 7. 

Equation (1) automatically insures an asymptotic approach of r 
to its optimal value 79, regardless of initial conditions. That the 
approach is always asymptotic can be seen from the following 
consideration. 

Let m be the lowest power in the development of F(z) into a 
power series at = 0. Since F(z) is monotonically increasing and 
F(0) = 0, therefore m is an odd number and the term is positive. 
Let n be the lowest non-zero power in the development of C around 
the point r= 79, for which C has a maximum, and hence for which 
dC/dr=0. Then n is even, and the sign of the term is negative so 
that C = C(ro) - g? (r-:79)". The development of dC/dr around 
7™ begins therefore with the term of the form —ngq" (r —7o)"~), 


C 
Therefore, the lowest power in developing F;(r) = r(%) into a 
7” |} 


power series of r at the point r = 79 is an odd number 2k + 1. Hence, 
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if A is a constant, in the neighborhood of r = rp equation (1) becomes 

dr 

dt 

Integrating this and denoting by r, the initial value of r for ¢ = 0 
we find 


=A Ces (2) 


(7; — 70)?* ae 
142A k(t, —179)2% ¢] - ’ 


which shows that r tends to ro asymptotically. 

It is, of course, possible to assume that accommodation does 
not maximize anything, but rather results in an adjustment to a 
particular value of something which varies monotonically with r. 
Effects of chromatic and spherical aberrations are, in general, 
such that for the same amount of ‘‘over-focussing’’ or ‘‘under- 
focussing’’ the detailed structure of the blurred image on the retina 
is different. The effect of chromatic aberration would seem to be 
ruled out by the circumstance that we can accommodate with mono- 
chromatic light in darkness. Effects of other aberrations cannot be 
excluded. We shall not consider them here, but their possibility 
should be kept in mind. It must, however, be remarked that the re- 
sponse to the structure of the blur would involve a much more 
complicated mechanism, requiring a response to a continuously 
varying ‘‘Gestalt.’’ 

A neurobiophysical interpretation of equation (1) presents, how- 
ever, some difficulties. As in loc. cit. (1951), p. 35, we may 
imagine that the contrast acts as a stimulus upon a pathway / 
which has the following characteristics: 


r=fo + 


(3) 


lepeatle Reon laarer (4) 


Q |p 


B 
a 


Then, every time there is a sufficiently sudden increase in C 
the pathway acts for a short time as an excitatory element, whereas 
for sufficiently sudden decreases in C, it acts as an inhibitory 
element. If now 7 is below the optimal value 7o and increases in 
sudden steps, the result is a succession of excitatory impulses in 
the pathway /, which, then through a proper connection, will en- 
hance the further increase of r. On the other hand, if r> 7 and 
increases, the C will decrease and will result in an inhibitory ac- 
tion of /, thus inhibiting the increase of 7, But if r<7o and ac- 
cidentally decreases, this will result again in an inhibition, and if 
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the inhibition 4nhibits the increase of 7, it will facilitate its de- 
crease, and we obtain an unstable condition. Hence this picture 
which is basically the same as that described on page 35 of loc. 
cit. (1951) must be further complicated. The complication arises 
from the need to have an adjustment to a maximum value rather 
than to any particular absolute value of C. On each side of the 
maximum, there are always two equal values of C, to which there 
correspond two different values of r. For the same value of C, r 
must either increase or decrease, depending on which ‘‘side’’ of 
the maximum that value of C lies. 


FIGURE 1 


Inasmuch as the mechanism of accommodation involves a large 
number of muscle fibers of the ciliary muscle and, therefore, also 
involves a large number of nerve fibers and neurones, we shall 
discuss the problem in terms of neuroelements and pathways, 
rather than in terms of individual neurones. One description can 
readily be translated into the other (Rashevsky, 1948, chaps. xxx 
and xlvi). 

Let us consider the following mechanism. 

First of all we do assume again that an increase of C results in 
an excitation of some neuroelement, whereas a decrease of C re- 
sults in an inhibition of the same neuroelement. This can again be 
achieved by introducing a neuroelement / the parameters of which 
are characterized by (4). In this case both excitation and inhibi- 
tion are produced by the same neuroelement J. By means of a cir- 
cuit shown in Figure 1, we can separate the excitation and inhibi- 
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tion spatially. Let / divide in two branches, /’.and /’’. The branch 
I connects at s, to a chain of excitatory neuroelements A1, Agere, 
which end with an excitatory link e. Let A and B be permanently 
excited centers consisting of self-circuited elements. Let A ex- 
cite a chain 6,, bg... of excitatory elements which end with an 
inhibitory link z. Let B, on the other hand, excite an excitation 
element C which connects at s3 to an inhibitory element d. The 
latter inhibits the chain 6;, bo,... at so. However, the branch /’ 
also connects at s3. 

When the contrast C increases, a short excitation occurs in /, 
I’, and /’’, Through s;, this excitation travels down the chain ay, 
@2,... to the element e, which becomes excited. 

At s3, however, the excitation only enhances that of d, and thus 
increases the inhibition at sg. Hence the chain b,, bo9,... and, in 
particular, its last link z remain unexcited. 

When C decreases, there is a short inhibition in J, I’, and /’’, 
The connection s; being inhibited, the chain a1, a2,... is not ex- 
cited, and e is also not excited. But the short inhibition of /’’ in- 
hibits s3, and therefore d ceases to be excited for a short time. 
As a result the chain 0,, b2,... becomes excited for a short time, 
and z, the inhibitory end link, also becomes excited for a short time. 

Thus now we have either a short excitation produced by e or a 
short inhibition produced by z (but never both simultaneously), de- 
pending on whether C increases or decreases. 

Now consider the following neural mechanism. Let L (Figure 2) 
be a center consisting of a large number of self-circuited neurons 
with thresholds distributed over a wide range of values. The total 
number of neuronic circuits which are excited at a given time de- 
termines the total excitation of L. Let that total excitation con- 
trol the contraction of the ciliary muscle. Then the greater the ex- 
citation, the greater the contraction of the ciliary muscle, and 
therefore, the greater the thickness 7 of the lens. 

Let 1 and 2 be two neuroelements which may be excited with in- 
tensities ¢«, and eg either separately or simultaneously. Let the 
excitation of 1 lead through 1’ to the excitation of 1°’, which acts 
on L. When 1’ is excited, the number of excited neuronic circuits 
in L will increase, and for a constant excitation of 1°’ the increase 
will be approximately linear with time (Rashevsky, 1948, chap. 
xxx, pp. 363 ff). Thus excitation of 1 results in an increase of 
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a-CHAIN b-CHAIN 


FIGURE 2 


excitation of L and therefore in an increase of r. The rate of in- 
crease dr/dt is approximately a linear function of «,, 

Let the excitation of 2 lead through 2’ to the excitation of the 
inhibitory pathway 2’’, which acts on L, inhibiting it. During this 
inhibition, the number of neuronic circuits in L will decrease, and 
hence r will decrease, the rate of decrease being approximately a 
linear function of €9. 

Neglecting thresholds we may consider the rate of increase of r 


as proportional to «;, the rate of decrease as proportional to eo. 
Putting therefore 


b = 17g (5) 
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and denoting by p” a positive constant we have 


dr i 
pre go (6) 

Let the chains 1, 1’, and 2, 2’’, cross inhibit each other as 
shown in Figure 2, and let the parameters be so chosen that 1’ and 
2° are not excited when e; ~¢g = 6 = 0. Then the difference of ex- 
citations of 1” and 2’, or of 1” and 2” is still proportional to 4, 
and equation (6) holds. 

Let the a and 6 chains in Figure 1 each branch in two, and let 
one element e and one element z connect with 1 (Figure 2), while 
another element e and another z connect with 2. Let a branch 3 of 
1’ inhibit one branch of the @ chain and one branch of the 6 chain, 
at Sg and sg correspondingly. Let a similar branch 4 of 2’. inhibit 
the other branches of the @ and 6 chains at ss; and s;. Since 1’-is 
excited only when ¢@ > 0, while 2’ is excited only if d < 0, there- 
fore when ¢ > 0, there are no effects of either e or z on the chain 
2, 2°, 2°. On the other hand, ¢ < 0, there is no effect of either e 
or z on the chain 1, 1°,-1°. 

Let us denote the absence of excitation of the elements e in 
Figure 1 by @; the absence of excitation of the inhibitory element 
~ in Figure 1 by 7. We now have the following situation. 

Casel. r> 70: 


dr dC : 
Sa) eat 2 Oe ees vs 
? dt dt : ? @) 
dr dC = 
24 NI 5 ees ee 5 | Disa Vee aan Sy § amit asia fay 8 
s dt at (8) 
Case II. r< 7: 
foe we &E 9 
> — Er ;. 
ee Gee aa es a; (9) 
d dC 
Pewee ake tt) ae e052 Bae, (10) 
dt dt 


In the case of equation (7), only z is present, but it does not af- 
fect the chain 2, 2’, 2’’ because of inhibition at s,. But s7 is not 
inhibited, and therefore there is additional inhibition at 1, de- 
creasing ¢. This will result in a decrease of dr/dt. With proper 
inequalities, which shall be discussed below, holding between the 
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parameters, the inhibition produced by 7 at 1 will change the sign 
of #, and the situation described by (8) will come into being. 
This means that now e will be present but not 7. Moreover, since 
for @ < 0 the connections s; and s7 are inhibited, but the connec- 
tion sg is not, d will tend to become more negative, resulting in a 
more rapid decrease of 7. Hence in case I, when r> 19, if r acci- 
dentally begins to increase, this increase will be immediately re- 
versed. If it decreases, the decrease will be enhanced. 

In case II, when ¢ > 0, an excitation e results, according to (9). 
Moreover for ¢ > 0 the connection sg is inhibited and e has no ef- 
fect on 2, but it does enhance center 1. Hence ¢ becomes more 
positive and a further increase of 7 results. 

If in case II d < 0, then z is present, according to (10). But this 
time, s7 is inhibited, and therefore z affects only 2, but not 1. 
Hence ¢ will become less negative. Under conditions mentioned 
above, and discussed below, ¢ actually will become positive, thus 
changing the situation to that described by equation (9). 

Hence if r< ro, any accidental decrease of r will be stopped, 
and reversed, whereas any increase will be enhanced. 

Thus in all cases r tends to r9, which represents the only stable 
configuration. 

An essential assumption implied in the discussion of the present 
model is that any accidental changes of « or j occur sufficiently 
suddenly, for the following reason. A slow variation of ¢ or 7 will 
result in a slow variation of 7, and hence of C, and such slow vari- 
ations will not result in temporary excitation or temporary inhibi- 
tion in the mechanism of Figure 1. 

While this assumption seems rather artificial, it is not altogether 
implausible. The mechanism has been described here in terms of 
neuroelements and pathways, each consisting of a large number of 
individual neurones or fibers. The continuity of the equations 
which govern. the behavior of neuroelements is only an apparent 
one. Actually the all or none law for each individual neurone 
makes all the continuous-looking changes step-like. There is, 
however, a problem here which will bear further elaboration. We 
could, of course, describe the whole mechanism in terms of indi- 
vidual neurones and fibers instead of a mixed neuroelement which 
is characterized by (4), and postulate the presence of on-fibers 
and off-fibers, which react to changes of C. 

We can, however, make a very rough evaluation of the conditions 
necessary in order that a small accidental positive change A@ in 
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case r> ro, or a small negative change in case r< rg would not 
only decrease, but actually result in an opposite change. 

Consider first case I, and let 6 = Adé>0. Then dr/dt = p? Ad 
according to (6). Since for r> 1, dC/dr< 0, this results, there- 
fore, in a decrease of C, the rate of decrease being given by 


dC dC ar dC a 
Paro a”: 4) 
This decrease of C will result, as we have seen, in a negative 
component A’ of ¢. Very roughly we may assume that A’¢é is 
proportional to dC/dt, in other words, because of (11) 
A’ 2 dC 2 2 dC A 12 
pete? eee AG (12) 
where g® is another positive constant. 
In order that the rate of change of r should reverse its sign, and 


; dC 
dr/dt become negative, the absolute value | p?q? = must be 
r 
greater than 1; or 
dC i. (13 
= > iS 
dr 2 q? ) 


The same result is obtained by considering in case II a small 
accidental negative increment of ¢. 

Inequality (13) cannot be satisfied for values of r which are 
too close to 79 because at r= 79, dC/dr=0. Hence the two roots, 
r’ and r’’ of the equation 


i: 
ga? 
give us the maximum obtainable accuracy of focussing, in terms of 
the constants p? and g?. For any value of 2, such that r° <r< 7’, 
the thickness r may fluctuate within that interval. As soon as, 
however, 7 becomes either less than 7” or greater than 7°” it is 
brought back into the interval. 

An electric circuit, ‘‘imitating’’ the hypothetical nerve mecha- 
nism discussed here, can, of course, be used for such devices as 
self-timing radio circuits, self-focussing cameras, etc. 


dC 


Se (14) 


The author is indebted to Dr. Herbert D. Landahl for a critical 
discussion of this paper. 
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Further consideration is given to assessing the degree of precision 
achieved in the localization of enzymes in cells by chemical staining 
techniques. Procedures discussed are those in which the enzyme acts 
on a substrate by a zero-order reaction, producing an intermediate which 
is then converted by an irreversible first-order reaction (velocity constant 
k) into a non-diffusible colored product. A relationship is obtained which 
permits the degree of localization for this type of process to be easily 
calculated. This localization factor is independent of the activity per 
unit volume of the site; it increases with the size of the site and, for ex- 
perimental time ranges, is virtually independent of the time of staining. 
It is shown that, under usual circumstances, satisfactory localization is 
achieved in sites of radius 1 p if k/D is of the order 109 cm™~? and that 
when this ratio has the value 1011 cm? localization is very good. 


In a previous paper (O’Sullivan, 1955), the general principles 
covering the localization of enzymes in tissues and cells were 
briefly discussed. It was shown that completely false localiza- 
tions could be produced in unsatisfactory methods and that careful 
examination of each method is necessary to ensure that real locali- 
zations are obtained. Skill in application of the methods is also 
of importance as, for example, inadequate preparation of the tissue 
slices can permit diffusion of the enzyme and produce erroneous 
results. The importance of the accurate mapping of cellular compo- 
nents and the possible introduction of quantitative applications of 
cytochemical staining make it highly desirable to assess, in addi- 
tion, the absolute precision of various possible methods. It has 
been suggested (O’Sullivan, Joc. cit.) that the ratio of the mass of 
dye deposited in the enzyme site to the total mass of dye produced 
by the action of the enzyme in the site during a given period would 
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provide the basis for a suitable measure of the accuracy of the 
process. For practical reasons, however, this was slightly modi- 
fied. The total mass of deposit has been replaced by the mass 
that would be produced in the given time if the intermediate, formed 
by the enzymic action, were completely converted into the deposited 
dye. This gave the localization index or factor F. 

Processes under consideration are those in which the prepared 
tissue is immersed in a solution containing a chromogenic sub- 
strate with a developing agent, both being present in excess. They 
rapidly diffuse into the tissue slice, and at sites of enzyme action 
the substrate is converted into a diffusible intermediate (diffusion 
coefficient D), which can react with the developing agent to give 
the dye deposit. The initial reaction is zero order of rate m 
(mol/cc/sec), and the final reaction is first order of velocity con- 
stant k&. The dye is taken to be non-diffusible. 

For the purposes of calculating F theoretically, a suitable 
enzyme distribution has to be selected, and the most satisfactory 
one would appear to be a uniform distribution in a sphere 0 &r< a, 
set in an infinite medium. For this system it has been shown 
(O’Sullivan, Joc. cit.) that the density distribution (mol/cc) of de- 
posited dye is given by 


2mnk * 
Cy = —! / $,(Du? + k)jt-1+ exp {-(Du? + Belldu, (1) 
0 


aT 


where 
Sin ru (sin au — au cos au) 
=> 
: u*(Du? +k)? 4 


and n, molecules of dye are produced from each molecule of dif- 
fusing substance. The localization factor is given by 


An [ ro jdr 


i= . 2) 


4 3 
rede mnt 


The further result previously given for F (O’Sullivan, loc. cit., 
Pp. 254) is erroneous having been obtained by the unjustified re- 
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versal of integration operations. Equation (1) is more conveniently 
expressed as 
Qmn,ka’ es F 
= | Pol(v +r°)p" — 1 + exp {—(v? +r?)u? de, (3) 
0 


where 
sin pv (sin v— v cos v) 
V=au,p=t/a,* =ka?/D, and p” = Dt/a? . 


Considering the terms in square brackets in equation (3), it is 
clear that when the dimensionless product \°? = kt is large, the 
other terms may be neglected in comparison with (v? + A”)y?. The 
percentage error in this approximation is always less than 100 
(kt)~' and is usually considerably less than this figure. As, in any 
reasonable cytochemical method, #¢> 100, this approximation is 
valid. In other words, for usual periods of staining and when & is 
not too small, the production of dye from the diffusible precursor at 
steady-state concentration predominates over the transient state 


production. 
Thus 
Qmn ,ka* t ; 
wes arD : 
where 
© sin pv (sin v — v cos v) dv 
ie J v?(v? + dr?) 


0 
can be evaluated by considering the function 


(1 - 2) {exp (1 — p) 22 — exp (1 + p) ait 
927(2? +r) 


in the 2(= 2+ ¢y)-plane. It is found that the integral of this func- 
tion from 2 = —« to = +~ along the z-axis is equal to 77 times the 
residue at 2=0 plus 272 times the residue at 2=%. Thus I 
becomes 


7 


as fpr -(A+ 1) exp (—A)-sinh pr} 
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and consequently 
mn ,t 


gw af on =(Nit 1) exp (—2\) -einkap Ate (4) 
d pr 


From equations (2) and (4) we get 
F~ 1-872 (A +1) exp (~ A)-(A cosh A - sinh 4). (5) 


This expression for F is the most convenient for use when A <1. 
When A is very small it can be seen from equation (5) that F ~ Nes 
The following form is more suitable for calculations when A > 1: 


Feri- = 17 Aya Ae 1) expe oa 

The localization factor F has value unity for complete localization 
and a value approaching zero for very poor localizations. It is seen 
from equation (5) that F is independent of time provided the in- 
cubation time is of reasonable duration. Prolonged staining times, 
however, will result in an intense deposit of dye which will signifi- 
cantly affect access of substrate to the enzyme and, clearly, com- 
puted values for F will be unreliable under these conditions. The 
factor is also independent of the specific activity (activity per unit 
volume) of the site. It increases with increasing values of & and a 
and with diminishing values of D. General experience with diffu- 
sion problems of this type indicates also that F is not greatly 
modified for sites of the same volume and alternative shape, pro- 
vided the surface area does not change considerably. 


TABLE I 


r 1000 100 50 10 5 1 On 0.01 
F 0.999 0.985 0.970 0.851 0.712 0.190 0.014 1074 
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Table I shows the variation of F with the dimensionless quantity 
\ and table II shows the way F varies with & (sec”') and D (cm? 
sec”') for a= 107* cm (corresponding roughly to the limit of reso- 
lution of the light microscope). Diffusion coefficients of artificial 
organic substrates are likely to be in the range given in Table II. 
Thus, if good localization is to be assured in small sites, the ratio 
k/D should possess a value near 10‘? cm~?, although &/D = 
10° cm? gives adequate results. Table III illustrates the way the 
localization factor varies with the radius of the site. The values 
given are calculated for D = 107° and with the & values given in 
the Table. 


TABLE Ii 


a i0-? ORES 105° 
aac ek) 0.190 0.851 0.985 
F (k = 10?) 0.572 0.953 0.995 
Some processes probably give poor localization due to the solu- 
bility and diffusibility of the dye which may result, for example, in 
the dye crystallizing in lipoid droplets. Mathematical assessment of 
the error arising from this source is also possible. 


Finally it is possible that the results quoted here have other in- 
teresting biological applications. 


The author would like to thank Dr. S. J. Holt for many discus- 
sions on cytochemical problems. 
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The problem of exchange of substances between circulating blood and 
tissue is discussed in terms of a model with continuously distributed 
tissue characteristics. Consideration is given to the distribution of 
perfusion rates, the distribution of metabolism coefficients for a sub- 
stance consumed at a rate proportional to its concentration, and the 
distribution of permeability coefficients for a slowly penetrating inert 
substance. Equations are derived for obtaining crude estimates of the 
moments of these distributions from exchange data. 


Many problems in applied physiology require a knowledge of the 
rate at which specific substances distribute themselves between 
the tissues and the blood. This knowledge is essential, for ex- 
ample, to a detailed understanding of the mechanism of decom- 
pression sickness, depth of anesthesia, and gas metabolism. On 
the other hand, this rate is frequently measured where concern for 
the specific substance is only incidental because the gas serves 
primarily as an indicator of some properties of the blood circula- 
tion. It is clear that the underlying theoretical framework is es- 
sentially the same for both of these interests. 

To cope with the problem of inert gas exchange, Behnke (1945) 
has fit his data on nitrogen elimination by the sum of two exponen- 
tial curves. He interprets this in terms of a model that lumps the 
body into a fat and lean body component. Using radioactive gases 
Jones (1950) found the sum of two exponentials to be inadequate 
and, stressing the fact that perfusion rates vary considerably from 
tissue to tissue, has used the sum of five exponentials to de- 
scribe his model of five lumped body components. 

Implicit in the above considerations are the assumptions that all 
the tissues are in parallel, that diffusion takes place very rapidly 
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in comparison to blood flow rates, (i.e., the entire process is 
limited only by the rate of blood perfusion), that the arterial con- 
centration is constant throughout the experiment, and that circula- 
tion times are negligible. 

A much more extended analysis of'the problem of gas exchange 
has been published by Smith and Morales (1944a,b) and Morales 
and Smith (1944; 1945a,b). In their theory, these authors have ac- 
counted for the effects of permeability, circulatory rates, and cir- 
culatory times. They also discuss the geometrical arrangement of 
tissues, but, for the most part, consider the tissues of the entire 
body to be in what they term a ‘‘competitive parallel arrangement.” 
(If diffusion and permeability are unimportant, the distinction 
vanishes between a ‘‘competitive parallel arrangement’’ and what 
is ordinarily thought of simply as a parallel arrangement.) Their 
results for the exchange rate is derived for an unspecified finite 
number of tissues and hence their general expression is in terms 
of an unspecified number of exponential decays. The decay con- 
stants are relatively complicated functions of the tissue param- 
eters. In applying the theory to the whole body, Morales and 
Smith (1945b) have simplified their expressions to the case of two 
types of tissues arranged in parallel with the same amount of blood 
flowing through each component per unit time. 

In contrast to Morales and Smith, Jones (loc. cit.) adopts the 
view that the effects of permeability and diffusion are negligible 
when compared to the effect of blood flow rate, but that the flow 
rate varies considerably from tissue to tissue. The variation in 
flow rate is borne out by comparing the actual values for various 
organs. To justify his assumption on the relative unimportance of 
diffusion or permeability, at least for the inert gases, Jones cites 
his observation that the difference in exchange rates of several 
inert gases with widely different diffusion coefficients can be ac- 
counted for solely in terms of their solubilities. This view has 
found further support in the theoretical deductions of Roughton 
(1952). 

In this paper we shall emphasize the variability of perfusion 
rates. That is, if we divide a gross tissue into many regions, each 
being supplied by its own characteristic capillary bed, then each 
region will have its own perfusion rate (i.e., flow rate per unit 
volume of tissue), and there is no reason to suppose all these 
perfusion rates to be the same. In general, there will be a distri- 
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bution of perfusion rates such that for a given rate there will be a 
certain number of regions that are perfused at that rate. These re« 
gions will have a certain total volume, and thus the volume of tis- 
sue that has a given perfusion rate is a defined function of that 
rate. We shall relate this function to the blood-tissue exchange 
rate with the hope of deducing properties of the function from ex- 
change data. 

The problem may be rephrased in slightly different terms. Thus, 
Supposing we randomly select a small volume of tissue and ask: 


What is the probability that that volume has a given perfusion 


rate? This probability will be given by the total volume of small 
regions that have that perfusion rate, divided by the whole volume 
of the gross organ. Hence the problem is to find a probability 
density. In some cases it will be seen how, in principle, explicit 
expressions which approximate this probability density may be ob- 
tained from blood-tissue exchange data. However in most in- 


stances we shall content ourselves with a determination of the 
parameters of this density (e.g., the mean, variance, and higher 


moments). We shall illustrate the method by the simplest case and 
then proceed to apply it to more complicated situations. 


Case I. Rapidly Penetrating Inert Substance. Single Homoge- 


-~ nous Tissue. Consider a substance X whose concentration in the 


arterial blood of a given gross tissue is suddenly changed from 
Cy to a new constant value C. Subdivide the tissue into N regions 
and consider the case where the solubility of X is the same in all 
these regions. Let 2, be the concentration of X in the zth tissue 
region, and y; its concentration in the blood leaving that region. 
If 7; is the total rate of blood flow through the region, then r,C’ 
equals the rate of flow of substance to the region, and r;y; equals 
its rate of flow from the region. The difference between these 
rates must equal the rate at which X accumulates in the region. 
That is, if V; is the volume of the region, 

V dx; 

Dae 

If diffusion and permeability effects are negligible, then the X in 
the blood leaving the region must be in equilibrium with the tissue 
X. Denoting the partition coefficient between blood and tissue by 


= 7;C —PiVi (1) 


a, this means that 
; Yi, = AX Ljep (2) 
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Substituting (2) into (1) and introducing the perfusion rate Aj; (i.e., 
A; = 7,/V,) we have 

az; 


=A,C — AjHay. (3) 
dt 


Considering 2; to be a function of both ¢ and Aj, the solution of 
(3) with the initial condition 2; = Co/& at ¢ = 0 is given by 


C 1 on 
#i(Aist) = — + 5 (Co ~ Cle aXe. (4) 


The amount of X present in the zth region at the start of the ex- 
C : 3 

periment is given by V; ae , and the amount present at time ¢ 1s 
Oo 


given by V; 2x; (Ai, t). The difference between these two quan- 
tities represents the contribution of the ith region to the total 
blood-tissue exchange: The total exchange is the sum of all these 
contributions. Hence, letting ¢(¢) represent the amount of X taken 
up by the tissue from the environment, we have 


N 


N 
b(t) = aE Vi Ere t) - “] eee ye V3 (2 eo ae) 
= 


a a 
int 

Since N is very large, we shall pass to the continuous case in 
the usual manner by taking the limit of (5) as N gets larger and 
each subdivision smaller. The order in which the summation is 
carried out is immaterial and hence we may choose first to sum 
over all tissues with the same X,, and then over all A;. Letting 
V(A)dA be the sum of the volumes of all tissue regions with per- 
fusion rate between A and A + dA, we may rewrite (5) in the form 


Cmte Pee 
$(t) = eee V(A) (1 - e72™*) dd 
0 (6) 
VolC-Co] C-Co £* 
aie ide eee Seal V(aye7 Man, 
a a 
0 
: (C - Co) 
where Vo is the volume of the gross tissue. Note that Vo ————— 
O 


is equal to d(), the asymptotic value of the tissue uptake, and 
can be determined directly by the exchange data. Note also that 


a ee Re 
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the last integral in (6) is simply the Laplace transform L(a 2) of 
the distribution function V(d). Rearranging (6) into the form 


a 


La2) =18~)- sl, (7) 


it becomes apparent that the deviation of the tissue uptake from 
its asymptotic value is proportional to the Laplace transform of 
the unknown distribution. Under practically all conditions of in- 
terest, expression (7) can be inverted to give an explicit expres- 
sion for V(A). To do this entails passing a continuous empirical 
curve, which can easily be inverted, through the measured points 
of d() — d(¢). Since the numerical computations in such a pro- 
cedure are rather laborious, and in anticipation of more compli- 
cated situations, we shall turn our attention to a slightly different 
approach. 

From (6), the nth derivative of d with respect to ¢ is given at 
t =.0 by 


co 


dd — : 
(4) -e (C co) | AV(A)AA. (8) 


The integral in (8) is simply the nth moment of V(A) about the ori- 
gin. Denoting that moment by p,, we have 


n+l oa 
mseqie 


gS te 9 
En = a“ 1(C — Co) ’ (9) 


and thus all the moments about the origin can be determined. In 
particular, the first moment »,, the mean of the distribution, is 
equal to the total functional rate of circulation through the whole 
tissue. 

The variance vp and higher moments v, about the mean can 
always be found from the p,’s by the relation (Mood, 1950) 


m= 2 (F) Candee’ (10) 


i=0 


If instead of cumulative tissue uptake our data are given in 
terms of venous concentration y, we simply replace (6) by the 
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expression 


ip OAa(A t)V(A)dr 
je (11) 


i “AV (aad 


0 


and proceed in the same manner as above. 

Thus from the exchange data alone, we can compute all the 
parameters of our unknown distribution. Under fairly general con- 
ditions, these parameters uniquely determine the distribution. For 
many purposes, the explicit form of the distribution is nonessential 
because information about the qualitative features (e.g., mean 
value, variability, and skewness) can be obtained directly from 
the parameters. However it must be borne in mind, in both this 
and in the following discussions, that the estimates of the param- 
eters are only as good as the assumptions involved in their der- 
ivation and, further, that even if the assumptions are very accu- 
rate, the accuracy of our estimates will depend on the accuracy 
to which the derivatives of the exchange data can be measured. 
Instead of determining the values of the derivatives directly from 
the data at time ¢ = 0, it is suggested that any convenient enipiri- 
cal equation be fitted to the exchange data and the required deriva- 
tives be calculated from this equation. In this way the determina- 
tion of higher derivatives from graphs are avoided, and the calcu- 
lated values will depend on the entire exchange curve instead of 
wholly on the initial process. 


Let us summarize the assumptions involved in the above ex- 
ample. First note that the use of the same constant C for arterial 
concentrations in all tissue regions implies that the tissues are in 
parallel. Equation (2) requires that the blood leaving each tissue 
region be in equilibrium with the tissue, or, in other words, that 
the exchange process is limited only by circulatory rates. The ab- 
sence of a subscript on the 4 in equation (2) implies that the 
solubility of the exchanging substance is uniform throughout the 
tissue. Also, the small durations required for the blood to travel 
from the arterial site of injection to the venous site of analysis 
have been neglected. 
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Case Il. Rapidly Penetrating Inert Substance. Heterogenous 
Tissue. The above results can be generalized to the case of non- 
uniform solubilities and a time dependent arterial concentration 
C(t) with little difficulty. Thus consider the case of M parallel 
regions, each with its own characteristic solubility, Oe ed ane, 
and each with its own characteristic distribution of perfusion 
rates V (A), j= 1,...,M. If the nth moment about the origin of the 
jth distribution is denoted by p,,, then by simply following the 
above outlined method, we can arrive at expressions corresponding 
to equations (8) and (11). 

Thus in place of (8), we now have 


and M if 
= Me aah {¢ oj)" *[C(0) --Colpnj + 


j=1 
ees ea a 
(-ayit ( =i) Biz sp (12) 
me ui) ) 


and instead of (11), 
ie 
dt"/ 4 
M we : d™-*C 
=e {-aprteo — Colun+1,j + a Gass) rer } 


fot 
M 
yes 
j=1 


As an illustration of the use of these expressions, consider the 
case of M = 2 (i.e., two types of tissue: fat and lean.) Expression 


(12) for M = 2, n = 1 becomes 


(=) Sy COy=C0)s (14) 
dt /+-0 


(13) 


and for n = 2, 


2 ac 
‘o = ~(Oipor + Oop122)LC(0) — Col + (ui1 + mor(G) (15) 


istributi i tes is given by 
The total volume distribution of perfusion ra 
V =V,+Vo. It is clear by definition that if uy is the Nth moment 
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of V about the origin then 

HEN =PN1 + EN; (16) 
and hence from (14) the mean p; of V(A) is given by 


dp 
(a). (17) 


Hi = fit + P12 = C(0) — Co" 


If two distinct substances with different solubilities are meas- 
ured simultaneously, then (15) will provide two simultaneous ex- 
pressions, each containing the same moments but with different 
o«’s and different concentrations. These two equations can be 
solved for the two unknowns po; and poe. Then using (14), (10) 
and (16) we arrive at the following expression for the variance of 
the V distribution: 


ve = fe + (uo —2)ur” (18) 
= (wor + 22) + (uo — 2) (wis + H12)"> 


where po = Vo the total volume of the gross tissue. Moments of 
higher order may be found by the same methods. 


Case Ill. Rapidly Penetrating Metabolized Substance. Let us 
turn our attention to a substance X which is consumed by the tis- 
sues. As before, we shall assume that X is held in simple solu- 
tion by the blood and tissues, and hence our results do not apply 
to such substances as O, and CO,. Consider the case where X is 
consumed at a rate proportional to its concentration z.. The pro- 
portionality constant per unit volume &; will be allowed to vary 
from tissue to tissue and from region to region within the same tis- 
sue. For example, &; may be a measure of the concentration of 
some enzyme involved in the consumption of X. In the ith tissue 
region, equation (3) will be replaced by 


ax; 
dt 


=AC -—A;Aa; - kya; (19) 


For simplicity, consider the case where the arterial concentra- 
tion is suddenly changed from Cp to the constant value C at ¢=0. 


Defining W (A,k) dAdk to be that volume of tissue which has a per- 
fusion rate between A and A + dA and a consumption rate constant 
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between & and k + dk, an expression for (é) may be written as 
fos) [o) t 
b(t) = J J {ext Berk Oyae E aka sR) (20) 


Differentiating the above expression n + 1 times with respect to 
t, and using relation (19), we find 


d™*'¢ 
oe = ap [ 


By solving (19) and differentiating the solution n times, d”x/ dt” 
may be found. If this expression is substituted into (21), for the 
time ¢ = 0, the result will be: 


= 
dgnt} = 


(=a)y>** > ae {|(")es ate Je] Nee eS em P (22) 


m=0 0 0 


(21) 


cones} W(A,k)dddk. 


If the probability of finding a tissue with a given k& is independent 
of its perfusion rate, (22) may be simplified by the following argu- 
ment. Let w(A)dad be the probability of finding a region with perfu- 
sion rate between A and A:.+dA, and K(x) the probability of finding 
a region with a rate constant between & and k+ dk. Then, the 
probability of both events occurring simultaneously will be given 
by the product of the two independent probabilities wad and Kdk, or 


W(A,k)dAdk = Vow(A)K(k)ddk. (23) 


Inserting (23) into (22), we see that the variables are separable. 
Hence defining S(n, m) by the relation 


s(n, m)=0"-™ (n)ee > fe ‘ o lett see Be 


we may rewrite (22) in the form 


n+1 
(—) =(-1)"*! ls S(n, m) R™K(k)dk + 
ae (25) 


Cop fe k"K(k)dk. 
0) 
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Denoting the nth moment about the origin of the probability dis- 
tribution K(&) by nz we have, finally 


a"*'¢ n= 
ant) = (Say ree ye S(n, m)n m = Copinn- (26) 
m=0 


Note from (24) that the coefficients S(n, m) can be determined 
as in the previous case by the simultaneous measurement of an 


inert, rapidly penetrating substance. Hence starting with n= 1, 
the mean and higher moments of the ‘‘metabolism distribution’’ 
can be found. 


Case IV. Slowly Penetrating Inert Substance. To treat the prob- 
lem of the influence of permeability, let us utilize the model of 
Smith and Morales (1944a). Thus, consider the volume of blood in 
functional contact with the ith tissue region to be given by B; and 
the volume of tissue in the zth region by 7. The volume V; of the 
‘tissue region’? as employed in previous cases would then be 
given by 


V,=B,+ Lx. (27) 


Letting h; be the permeability coefficient and S; the surface 
area of the blood vessels in functional contact with the tissues in 
the zth region, the equations of Smith and Morales would read 


dx; 
ihe rr AiSi(yi -:% @3), (28) 
dy i 
B; “He TEE Yi AS (ys Oy), (29) 


As a first approximation, we shall restrict ourselves to those 
cases where the perfusion rates are primarily determined by 6,, 
the number of open capillaries in the ith region. Then the total 
rate of flow through that area will be proportional to 0; and also to 


6 P the pressure drop across the capillaries. Since 6; will be pro- 
portional to B; we may write 


6P, (30) 


where A is a proportionality constant. 
If A; is regulated primarily by 6;, then S; will be approximately 
equal to 0; 27p1, where p and 7 are the average radius and length 
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of the capillaries respectively. On the other hand, B; is equal to 
0 zp 7l so that 


OB; 
es 
Denoting B;/V ; by «; and ASP by y, we see from (30) that A; = €iy. 


Let us define the total ‘‘concentration’’ p; in the heterogenous 
tissue region of volume V; by the expression 


Si= (31) 


Biyi+ Tie; 
V; : 


Then using relations (30), (31), and (32), we may transform (28) 
and (29) into 


Pi = (32) 


ap (Cc 33 
dt Sey ks — Yi), ( ) 

and 
WY (Cy, ae (1 ] 34 
pene! i et a (Le; + es)yi — Op (34) 


In passing to the continuous case, we shall assume that y and p 
are relatively invariant with respect to different tissue regions. 
This is in accord with our assumption that flow rates are regu- 
lated primarily by the number of open capillaries. Further, for 
simplicity we shall make the plausible assumption that B; << T; 
and hence 1-¢«,;=1-8;/(B;+T7;)~1. Since d; = ye,, (33) and 
(34) can be rewritten as 


ap; 
at 


d E 2h; ar; 
=y(C --y;) - f: a y= 93 (36) 
dt p y 


Then the tissue uptake ¢(¢) will be given by 


= AC < Yi) (35) 


g(t) = i il Gc) Os aan (37) 
(@) (0) 


where V(A, 4) dAdh is the volume of tissue having perfusion rate 
between A and \:+ dA, and permeability between A and A + dh, and 


Po is the initial value of p. Differentiating n times with respect 
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to ¢ yields 


ie eles VG, h) dadh. (38) 


In order to apply (38), consider the case where at the time ¢ = 0, 
p=y=0, and the arterial concentration is suddenly increased 
from zero to a.constant value C. Using these values together with 


n 
(35) and (36), we can then find i . By substituting this 
t=0 
value into (38) and assuming that the distribution of permeability 
coefficients is independent of the distribution of perfusion rates, 
the desired relations between the moments of the two distributions 
can be found. Thus, for example, if &€, is the nth moment about 
the origin of the permeability distribution we have 


- 
RE oe = pil, (39) 
ele 
a*¢ 
= =—y pil, (40) 
- t=0 
dd 2 &1 
Ty ee (41) 
t=0 p 


d‘¢ &2 
(Gr) --[s seat (y? + 4y? = + ty ue. (42) 


If the perfusion rate moments p, are measured independently by 
use of a rapidly penetrating substance, then (40) can be used to 


1 


obtain y. The value of — and =a can then be obtained from (41) 
p p* 


and (42). On the other hand if the uz, have not been measured, LA 


can be obtained from (39), y from (40) and oa from (41). 
p 


I would like to thank Dr. Clifford S. Patlak for reading and dis- 
cussing this paper. 
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A non-constant-sum continuous game analyzed in a previous paper by 
one of the authors (A. Rapoport, Bull. Math. Biophysics, 18, 15-30, 1956y 
is extended from two to NW players with special emphasis on V = 38. It is 
shown that the concept of stability becomes in this case one of §‘pairwise 
stability,’’ and depends on the so-called distribution matrix of rewards. 
The distribution matrix and the collusion structure jointly determine the 
end states of the game. Conditions which lead to the emergence of one, 
two, or no ‘‘parasites’’ are derived. An apparatus is described which 
provides a physical analogue of the game, making possible the isolation 
of behavioral variables under the prescribed conditions of the game. 


The two person non-constant-sum continuous game considered by 
A. Rapoport (1956) can be generalized to an N-person game as fol- 
lows. The satisfaction (or utility or pay-off) of the zth player is 
given by 

N 
S,=log [1+ )° a,,2,]-B,2;, @=1, 2,--+, N) (1) 

j= 
where a,, is the fraction of the effort x, of the jth player received 
by the ith player. Rapoport (op. cit.) has investigated the case 


where @,, = Qo. =D} Gy = = GP+I=15 By =Bo- 
For N players, Rapoport’s optimal lines become optimal hyper- 
surfaces which, in the absence of the ‘‘social matrix,’’ are hyper- 


planes given by 
N 


a,,=B;,A+ Ne 5; @;) (2) 


jemi 
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The solution of the system (2) in general (except in degenerate 
cases) defines a point (7,, %,.-., %,) in the N-space of the ef- 
forts v, This is an equilibrium point in the sense that the satis- 
faction of each player is maximized at that point with respect to 
that player’s own effort; i.e., if all the other efforts a, (j #2) re- 
main constant, no player is motivated to depart from this maximiz- 
ing value %,. This is the Nash-Rashevsky equilibrium. In general, 
the maximizing attempt of each player is such that he tries to bring 
the point (#,, %,---, 2 ) onto his own optimal hyperplane. Whether 
or not these attempts result in the point approaching the Nash- 
Rashevsky equilibrium depends on the stability of the equilibrium. 

However, whereas in the two-person case the equilibrium is either 
stable or unstable, in the N-person case, stability must be examined 
with regard to each pair of players separately. Let the efforts of 
all players but two be fixed and consider the resulting two-person 
game between Players 1 and 2, namely, 


S, = log (1+ @,, 2%, + @,,%, + &,)-B, 2, , (3) 
S, = log (1+ @,, %, + Go %, + ky) — By % (4) 


where the &, are the lumped contributions to S, of the remaining 
players (considered fixed). The optimal lines, determined by 
58 ,/da, = 58,/8x, = 0 will be 


44, ~ By (1+ ky) — By (Gy, 1 + Oy %Q) = 0, (5) 
Ag ~ By (1 + kg) = Ba (@q, @ + Gy, 25) =0 « (6) 
We see that the slopes of the optimal lines in (%,, %)-Space are 
independent of both the efforts and the reluctance coefficients B ie 
The stability of a two-person game of this sort has been shown 
(Rapoport, op. cit.) to depend only on the relative absolute magni- 


tudes of the slopes of the optimal lines. Therefore, following this 
result, we have 


Theorem 1. Players ¢ and j are stable relative to each other if 
and only if 


a,,/0,,> 4,,/4,, OF @,,4,,— @,,4,,>0. (7) 


It follows that the pairwise stability of the total game depends en- 
tirely on the signs of the two-rowed principal minors of the distri- 
bution matrix (a,;)- 

The 3-person Game. In what follows, the players will be labeled 
X, Y, and Z, and their respective efforts, z, y, and zg. Let the frac- 


PARASITISM AND SYMBIOSIS 221 


tions of each effort distributed among the players be p, g, and 
t (p+ q+r=1), where it is understood that each player keeps a 
pth fraction of his own effort and distributes the remaining frac- 


tions among the other two. There are two types of distribution 
matrices: 


prg@ PQaq 
A,=|qprjand A, =|qpr|. (8) 
r gp rrp 


3 
Under the restrictions a,,= p and Bhs a;,= 1, every other matrix 
j=1 
reflects simply a re-naming of the players. We call the case repre- 
sented by A, the cyclic case and the other the non-cyclic. The 
term ‘‘cyclic’’ derives from the fact that the fractions qg and rf ro- 
tate cyclically among the players, that is, X gives a gth fraction of 
his effort to Y, who gives a qth fraction to Z, who gives a qth frac- 
tion to X, and similarly for 7. This is not the case in the non- 
cyclic case. Examining the principal minors of A,, we have 
immediately 
Theorem 2. If the distribution matrix of a 3-person game is 
cyclic, the game is either stable or unstable with regard to every 
patr of players. 
Note: this theorem does not hold for N > 3. 
Examining now the non-cyclic case in the light of Theorem 1, we 
see that the conditions for stability of the three pairs of players 
are as follows: 


X vs Y is stable if p —- g> 0; 
Y vs Z is stable if p — r> 0; 
Z vs X is stable if p? - gr>0. 


We see that stability (instability) of X vs. Y and of Y vs. Z implies 
the stability (instability) of Z vs. X. However, it is possible to 
have four distinct cases: (1) all pairs stable, (2) one pair stable, © 
(3) two pairs stable, (4) all pairs unstable. The conditions of 
stability are diagramatically illustrated in Figure 1. Each point 
within the triangles represents a possible pair of values (7; 2), 
which also determines r=1-—p-—g. The unshaded region is un- 
stable for all pairs; the vertically hatched regions are stable for 
one pair; the horizontally hatched regions are stable for two pairs; 
the cross hatched region is stable for all pairs. 
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(a) CYCLIC (b) NON-CYCLIC 


FIGURE 1. Stability in (p, 9g, ")-space with no collusions. Unshaded: 
no pair stable; vertical hatched: one pair stable; horizontal hatched: two 
pairs stable; double hatched: all pairs stable. 


We now deduce the end states of the game. If all pairs are 
stable, the Nash-Rashevsky equilibrium will be maintained. In the 
unstable cyclic case (all three pairs unstable), the only persistent 
end state is that of two ‘‘parasites’’ supported by one host. This 
is so, Since every two-person sub-game is unstable and results in 
one of the players becoming a parasite on the other (Rapoport, 
op. cit.). Who will be the host is a matter of chance. The non- 
cyclic case remains to be considered. Similar considerations lead 
to the following easily deducible results, referring to all possible 
situations (except for re-naming of players). 


1. If Y vs. Z and Z vs. X are stable, X vs. Y unstable, then 
either X is a parasite on Y and Z, or Y is a parasite on X and Z. 

2. If Y vs. Z is stable, X vs. Y and Z vs. X unstable, then 
either X is a parasite on Y and Z or Y and Z are parasites on X. 

3. If X vs. Y is stable, Y vs. Z and X vs. Z unstable, then 
either Z is a parasite on X and Y or X and Y are parasites on Z. 

4. If X vs. Y and X vs. Z are stable, and Y vs. Z unstable, then 
either Z is a parasite on X and Y, or Y is a parasite on X and Z. 


In summary, we have, without distinguishing the players, the 
following possible situations: 

1. All three ‘‘working’’ (stable case), 

2. One parasite (one unstable pair), 

3. Hither one or two parasites (two unstable pairs), 

4, Two parasites (all unstable pairs). 


Collusion. Collusion among a subset of players in an N-person 
game occurs when the members of the subset act essentially as one 
player, that is, choose a collective strategy. A generalized notion 
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of collusion was introduced by Rashevsky (1951) in the form of a 
matrix, which is essentially a linear transformation from the vector 
space of satisfactions to a vector space of **secondary’’ satisfac- 
tions. Thus the matrix introduces linkages among the players’ 
satisfactions. Rapoport has called it the ‘‘social matrix.”? The 
special case where the columns of the social matrix are all equal 
corresponds to the ‘‘grand collusion’’ of game theory, where all the 
players are in collusion against ‘‘nature.’’ All other collusions 
can be represented by special forms of the social matrix. But 
intermediate cases of ‘‘partial’’ collusion are also representable. 

Let the social matrix be (B,,)) that is, the matrix of the linear 
transformation 


1 
Spe 8,5) (9) 


Here Ss} is the ‘‘secondary”’ satisfaction of the zth player, and 
B,, measures the. dependence of S} on S,. The optimal hyper- 
surfaces (no longer hyperplanes) will then be given by 


yy i _f, = 0. @=1,.-+, 9) (10) 


Suppose now that these N hypersurfaces intersect at a point 
(Z,, Z,---, &,). To determine the stability of Players 1 and 2 
relative to each other, we consider 2, (j > 2) fixed so that equa- 
tions (10) for ¢=1, 2 define a pair of curves in (z,, 2, )-space, 
namely, 


F, (@,, @) = y ptren ~B,=0, (11) 
Pesca 9 Fe i oY 
N Beea- | 
3j 7? ~ B, = 0 (12) 
1 ORM (ree Or ed ie | 2 ’ 
2 (#1; 2) Lisa h ten Fh 


N 
ae " 
where k, = »s a,,v, are held constant. Stability, we have seen, 
i=3 
depends on the relative absolute magnitudes of the slopes at their 
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intersection. However, whereas in the non-collusion case these 
slopes were independent of the point of intersection (and therefore 
of the contributions of the other players), such may not be the case 
now. We compute the slopes by the implicit function method. 


(=) 5F,/8a, 
dx 1 oF, / 8x, a 
N 
ay By 51 /(1 + Gy, @y + yy @y + h,)* 
j=l 
= ’ 
N 
DZ By 5% 1 %o/(1 + Oj, By + Ayg Ze + k)° 
j=1 
a OF, /dx, 
day], OF, /da, hab 


N 
a By jg jy /(1 + 0,1 2, + Ojo ty + k,)? 


N 
= By ,%jg/(1 + ay, @, + G0, +k)? 


jut 


Therefore the condition of stability between Players 1 and 2 for a 
fixed contribution from the other players becomes 


N 

2 _ ee: 2 
We By U/L + a), By + Oj. B + hy) 
fat 


N 

ay By 5% Og /(1 + a,) 8, + a, Be + ke)? 

pat 

(15) 


= By {4% aj,/1+ 4,8, +43 % + k,)° 
j=l 


N 
Es ay /(1+a,,%, + Gio @ + k,)* 


ped 


PARASITISM AND SYMBIOSIS 225 


where z, and Z%, are solutions of (11) and (12). This inequality 
may or may not be satisfied for a given set of values of k, » which 
depend on the contributions of the other players. Hence it is not 
sufficient in this case to investigate the pairwise stability at some 
one point (as has been done in the case of the diagonal social 
matrix). Two players may be stable with respect to each other in 
some regions of the effort space and unstable in others. 

We have seen that the social matrix offers a generalization of 
collusion theory. In particular, when the entries are all either zero 
or unity and if the social matrix is a ‘‘direct sum’’ of two or more 
matrices with unit entries only, we have collusions as they have 
been described in game theory, i.e., coalescences of subsets of 
players into single players. 

We shall investigate the case of a complete collusion of two 
players against the third. The case subsumes the subcases of the 
cyclic and the non-cyclic distribution matrices, and in the latter 
case the various combinations, e.g. the collusion between a stable 
pair, an unstable pair, and the stability of the third player with re- 
spect to each of the players in collusion. 

While there exist three possible collusions of two players against 
the third, namely, X and Y vs. Z; X and Z vs. Y; and Y and Z vs. 
X, it is sufficient to study the first two for all possible sets of 
values of p, g, and 7, since the last collusion will then be reduced 
to one of the two investigated by a renaming of the players. More- 
over, in the cyclic case, the study of only one collusion exhausts 
all possible cases. In what follows, results which are self-evident 
are given without proofs. (See Figure 2.) 


I. Cyclic case: X andY vs. Z 
1. X vs. Y is stable if 


(p? — 1g)? > 0. (16) 


Hence the colluding pair is always stable (or neutral in the special | 


case where p? = rq). 
2. X vs. Z is stable if 


p(p? - gr) K2 + g(pq-) Ky > 9, (17) 


3 
Kg is pair is stable in the en- 
where K, = (1 + 2 a; @,) . Hence this p 


ail 
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tire octant w, y, 2 > Oif p’ > qr; pg > ?”, and unstable in the entire 
octant if p? < gn pq< r?, In the region of (p, g, f)-space where 
(p? — gr) (pq- r*) < 0, the stability of X vs. Z depends on the point 
(x, y, 2); i.e., there are regions in the octant 2, y, 2> 0, where the 
pair is stable and other regions where it is unstable. 


I. CYCLIC; X AND Y IN COLLUSION 


Gr 
Gt 
7, 


Xvs Y XvsZ 


ILa. NON-CYCLIC; X AND. Y IN COLLUSION 


IT b. NON-CYCLIC; X AND Z IN COLLUSION 


p p = p 
Xvs Y Xvs Z Yvs Z 


FIGURE 2. Stability in (p, 9g, r)-space with collusions. Unshaded: 
unconditional instability; hatched: conditional stability; shaded: uncon- 
ditional stability. 


3. Y vs. Z is stable if 


K? r(pr- gq’) + Ke p(p? - qr) > 0. (18) 


Hence this pair is unconditionally stable if pr — q? > 0, p? — gr> 0; 
unconditionally unstable if pr—q*® <0, p? ~gr<0; and condi- 
tionally stable otherwise. 
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Ila. Non-cyclic case: X and Y vs. Z. 

1. X vs. Y is unconditionally stable (or neutral) as in the cyclic 
case. 

2. X vs. Z is stable, unstable, or conditionally stable under the 
same conditions as in the cyclic case. 

3. Y vs. Z is stable if 


q (p-1) K? + p(p? - 17) K2 > 0. (19) 


Hence this pair is unconditionally stable or unstable, depending on 
whether p> rorp<r. 


Ilb. Non-cyclic case: X and Z vs. Y 
1. As before, the colluding pair X vs. Z is always stable. 
2. X vs. Y is stable if 


K? p(p? - g) + K3 1? (p- g)>0. (20) 


Hence this pair is unconditionally stable or unstable, depending on 
whether p> gorp< q. 
3. Y vs. Z is stable if 


K?2 @ (p—1)+K? p(p? - 17) > 0. (21) 


Hence this pair is unconditionally stable or unstable depending on 
whether p> rorp<r. 


Note the lack of symmetry between cases Ila and Ib. In the one 
case only one non-colluding pair is unconditionally stable or un- 
stable; in the other case both non-colluding pairs are uncondi- 
tionally stable or unstable. This is due to the genuine difference 
in the structures of the situations. The particular distribution 
matrix A, reveals a reciprocity between « and y (they exchange a 
qth fraction of their efforts). Hence a collusion of this pair leads 
to results different from those obtaining from the collusion of X 
and Z. The pair Y and Z, which has a reciprocity based upon the 
exchange of 7, shows a collusional stability. 


The Plane of Dependency. 
We have seen that in certain cases involving collusion, pairwise 


stability depends not only on the entries p, g, and r of the distribu- 
tion matrix but also on the efforts z, y, and z. It remains to de- 
termine this dependence. To fix ideas, take the case represented 
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by equation (17). Substituting for Kk, and K, and re-arranging, we 
find 


(p- gC) w+(q-pC)y+(q- 70) 2+1-CH=0, (22) 
where 
C =[p(p? - gpl”? (g(r? - pg)" ”?. (23) 


Equation (22) defines a plane in (x, y, 2)-space which separates 
the region of stability from that of instability. To have significance 
in our model, the expression (23) must be real, and the plane (22) 
must have an intercept in the octant a, y, x >0. The first condi- 
tion is satisfied already by the condition of dependency (i.e., that 
stability is conditional), for in that case (p* — gr) and (pg - r”) 
must have opposite signs and hence (p? = gr) and (r? — pg) must 
have like signs. The second condition will be satisfied if either 
one of the following inequalities hold: 


(p= gC) (tC) 0; 
(g- pe)(l-C)<0; (24) 
(gO) (0) 05 


“‘Sympar’? 

The problem of putting the foregoing theory to an experimental 
test involves a decision of just how much of the situation should 
be ‘‘structured.’’ One might, for example, build a servo-mechanism 
whose construction is a physical realization of the mathematical 
structure of the problem. Such a mechanism would necessarily be- 
have exactly as prescribed unless there are errors in the derivations 
or in the construction. Of course, the behavior of such a mechan- 
ism in no way constitutes a ‘‘test’’ of the theory except in the 
trivial (physical) sense. At the other extreme, situations might be 
contrived in which people or animals expend efforts to obtain re- 
wards which they share in accordance with self-imposed rules. It 
would be quixotic to expect that our equations would in any signif- 
icant measure represent this ‘‘realistic’? situation—there are too 
many arbitrary assumptions in the model. 

An intermediate strategy is to structure some of the situation 
and leave some to the ‘‘free will’? of experimental subjects. Ac- 
cordingly, there was constructed at the Mental Health Research 
Institute (University of Michigan) an apparatus labeled ‘‘Sympar’? 
(a contraction of ‘‘symbiosis and parasitism game’’). The distribu- 
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tion matrix and the form of the satisfaction functions S, are built 
into the machine (with arbitrary settings of the parameters p, q, 7, 
B,). Likewise an arbitrary social matrix B,, can be plugged in. 
However the efforts 2, y, 2, are left entirely to the subjects, who 
are only instructed to ‘‘try to maximize’’ the $,, which they ob- 
serve as pointer readings. The fact that the acts of each subject 
affect all of the S, and the subjects’ ignorance of the structure of 
the ‘‘game’”’ causes considerable confusion at least in the begin- 
ning. The object of the experiments is to determine whether and 
how the behavior of the subjects becomes stabilized as a result of 
experience, whether in the absence of communication among them 
the stabilization will occur near the Nash-Rashevsky equilibrium 
(if stable), whether in the unstable cases the ‘‘parasitic’’ relations 
predicted by the theory will be preserved, and especially whether 
when communication is permitted, the subjects are able to reach 
the ‘‘social optima’’ (whose stability depends on continued co- 
operation of the players and their ability to resist temptation to 
achieve greater gains at the expense of their partners.) 

‘*Sympar’’ is based upon the Sterling LM-10 Analogue Computer. 
(See Figure 3) 

A voltage representing the output (effort) of a player (X) is de- 
veloped across a potentiometer under the control of that player. 
Other potentiometers, controlled by the experimenter, take frac- 
tions of this output pz, gz, re and 8, a. By the use of a two posi- 
tion switch the experimenter selects the distribution matrix A, or A, 
that he wishes to investigate. The appropriate voltages are pre- 
sented to a summing amplifier which has an output V =1+ px+ 
gy +z. The logarithm of V is taken by use of a network of non- 
linear elements, that relies upon the current-voltage characteristics 
of germanium diodes. Due to the large losses entailed by the use 
of this network two stages of amplification are required before 
there is sufficient power to drive the output meters. Before the 
second stage of amplification the signals representing ‘‘log V”’ 
and **8,«’ are rectified and subtracted. The meter thus is pre- 
sented with a voltage proportional to ‘‘S, .”” 

It was quickly discovered that in order to increase the difference 
in reward between the Nash-Rashevsky equilibrium and the social 
optimum the constant ‘‘1’’ in the logarithm had to be increased. 
This was made adjustable (a in Figure 3) and its logarithm is sub- 
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tracted from the resultant function, thus giving instead of equation 


(1) 
a+ pet qyt ra 
s, = os é 2 )-B,2. (25) 
a 


In a two-person game the values p = 0.7, g= 0.3, B, = B, = 9.2, 
a=3 give almost a 2/1 ratio of social optimum to the Nash- 
Rashevsky equilibrium. This makes the reward for cooperation 
more ‘‘dramatic.’’ Initial trials of the apparatus, with the above 
constants, using students as subjects, tend to indicate that when 
verbal communication is permitted a rapid approach to the social 
optimum is made, but not otherwise. Indeed in the latter case, the 
efforts oscillate about, rather than approach, the Nash-Rashevsky 
point. Much more data must be collected and analyzed before 
credence can be given to this result. If it is confirmed, interesting 
possibilities are offered to control the level of group performance 
by varying the restrictions on its communication patterns. 
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An error in a previous paper (Bull. Math. Biophysics, 18, 31-56, 1956), 
pointed out by Professor R. L. Wilder, is corrected. 


In a previous paper entitled ‘‘The Geometrization of Biology”’ 
(Rashevsky, 1956), we discussed some special types of topologi- 
cal spaces which have the following characteristic: The space M 
is.a subspace of another space S. The loss of a point B of M 
makes some points 41, G2,..., 4; of S = M acquire with respect to 
M a property P which those points do not possess unless the point 
B is removed from M. The particular property P considered in loc. 
cit. is to make M connected. 

We have given in loc. cit. several examples of such spaces. 
Professor R. L. Wilder has pointed out to the author that one of 
the examples contains an error and is altogether wrong. Professor 
Wilder also showed the way to correct this error and to make the 
example valid, besides pointing out several minor errors. 

On page 44 of loc. cit. we considered a case of a space M, 
which we called a A*-space, and which is described as follows: 

Consider a family of circles in the plane E?: 


(a=)? +y? = 0, (1) 


where & = = or & is an irrational number between 0 and +. The 
cardinal number of such a set of circles is 8. Each circle inter- 
sects the line y =.0 at x =-0 and at x = 2a. Remove in each circle 
the point 2 =.2a. The space M obtained by this removal of points 
from the family of circles (1), considered as a subspace of 
S =-E?, is connected. Removal of the point 0, which we shall 


designate by 6, from M makes it not connected. 
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Following the above description we stated on page 44 of loc. 
cit. that ‘‘in order to make it connected with the point O removed, 
we must add to it all the points « =-20.’’ This is obviously wrong, 
for even with all points 2 =-20 present, the removal of the point O 
makes the family of circles homeomorph to a family of parallel 
straight lines, thus leaving it-non-connected. 

To make the example valid and to have the removal of the point 
QO make a non-countable number of points acquire property P, we 
proceed as follows: 

In the family of circles (1), let a =+ or & be equal to a rational 
number. The removal of the point B makes the family of circles 
(1) non-connected. But now the addition of all points of the seg- 
ment (0, 1) of the z-axis which do not lie on any circle of the fam- 
ily, that is, of all points with e420, 0<a#<1, y=-0; will make 
the resulting space connected. The number of points to be thus 
added is non-countable, being the number of all irrational points 
between 0 and 1. Thus the removal of the point B makes all 
points 2 4.20, 0<a<1, y =-0 acquire the property P. 

The A**-space described in paragraphs 4 and 5 on page 44 also 
does not satisfy the necessary conditions. Let z2=-2, be an ir- 
rational number, 0 < 2; <1. The portion of the space S defined, 
for which 2 > 2,, call A; the portion for which 2 < 2;, call B; and 
the portion for which z =-z;, with the point B of this set removed, 
call C. Then A is connected, being topologically similar to S. 
But the irrational number £ is a limit point of the irrationals B’- be- 
tween 6 and 1, and correspondingly every point of C is a limit 
point of A. Thus AUC is connected, as it consists of a con- 
nected set A plus limit points of A (Wilder, 1949, p. 16, Theorem 
7.2). Similarly B UC is connected. Finally, A UB U C is.con- 
nected, being the union of two connected sets with points in com- 
mon. Thus S-8 is not disconnected (/bid, Theorem 7.3). 

If the example is modified so that each \’-space (in the amended 
sense) is in a plane which passes through the z-axis, in a sequence 
of planes with each two adjacent planes forming an angle z/n, 
where n may be as large as we wish, then the point B =.(0, 0, 0) 
will be a disconnecting point, and the set of irrational points of 
the x-axis (0 < #< 1) will acquire property P. This space we may 
call \’*-space. 

The theorem stated on page 44 remains in force for those \’ and 
A spaces. However, for precision the theorem should read: ‘‘The 
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analytical intersection of two \(A’, \”%) spaces is not a A(A’, d’’) 
space and is always not connected.’’ 


The following minor errors should be corrected: 

Line 20, p. 42: The number of sets that can acquire property P 
is not countably infinite; it is the same as the number of ways of 
selecting subsets of the B;’s...i.e., 2%° =X. 

Line 5, p. 43: Analogous remark. 

Line 2, p. 44: The ‘‘set of sets’’ is of cardinal &., 

On this occasion we would like also to suggest a change in a 
sentence on page 47 of loc. cit., which, while not wrong, is am- 
biguous. The sentence on lines 7 to 10 from the top should read: 
‘*The removal of point 8; from M makes either the point , or the 
point &_ of S— M acquire the property P; while the removal of the 
point B2 makes either the point a! or the point a5 of S—M ac- 
quire the property P. If, however, after removal of 6,, a1 is 
added, then &» loses the property P and vice versa. The same 
hold about B2, 43, and a.” 

The author wishes to express his gratitude to Professor R. L. 
Wilder. 
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The idea of assigning an information content to a graph is extended to 
include the following two situations: (a) Combinations of sets of topolog- 
ically equivalent points are used as symbols; (b) Points of the graph may 
exist in different states. 


In a previous note (Trucco, 1956) we discussed the relation be- 
tween the group, H, of a graph, G, and the topological information 
content of the graph, as defined by N. Rashevsky (1955b). 

We consider graphs whose points may be of different categories 
or kinds. Each point belongs to one and only one category; as 
extreme cases the points can be assumed all alike or individually 
different. 

Two such graphs, G, and G,, are topologically congruent if there 
exists a one-to-one mapping of G, onto G, satisfying the condi- 
tions (1), (II), and (III) given on page 186 of G. Pélya’s paper 
(Pélya, 1937). In dealing with oriented graphs we may add to 
Pélya’s condition (I) that the mapping should preserve the orienta- 
tion of the lines. Congruence between graphs is a reflexive, sym- 
metrical, and transitive relation. 

All congruent self-mappings, i.e. all the one-to-one mappings of 
a graph G onto itself which satisfy those three conditions, form the 
group H of the graph. The transitive sets of H determine the sets 
(or classes) of topologically equivalent elements (points and lines) 


of G. 
237 
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From this definition [cf. condition (I1)] it follows in particular 
that two points of different kinds cannot be topologically equivalent.* 

The topological information content, or entropy per symbol, of 
the graph (due to its points) was defined by N. Rashevsky as 


I=-)° u, leu,» (1) 
k= 
where lg denotes the logarithm to the base 2, and 
n 
w,= —; (2) 
n 


m is the number of sets of topologically equivalent points in the 
graph, and n, the number of points belonging to the #th set. The 
graph has n points, so that 


> Mame (3) 


In general, we consider a discrete source of information (cf. 
Shannon and Weaver, 1949, pp. 9 ff.) as being described by a sim- 
ple Markov process: The source may exist in m different states, 
and a symbol or ‘‘letter’’ is emitted in the transition from state j to 
state &, which occurs with probability Pine 
We have 


m 
pm CORTE, 2m (4) 
k =i 


The states can be represented as points in a plane and any non- 
zero probability of transition, p,,, as a directed line going from 
point 7 to point k. In this way we obtain the source-graph, not to 
be confused with the other graphs dealt with here. 


*It should be mentioned that Pélya considers the graph as a chemical 
formula and sets up the group H for what he calls a chemical ‘*‘Stammkorper.”? 
This is essentially the same as A. C. Lunn and J. K. Senior’s (1929) 
chemical skeleton in which all monovalent substituents are supposed to 
be of the same kind (H-atoms in Pélya’s example). Thereupon different 
monovalent atoms or radicals are substituted onto the skeleton, but the 
number of possible isomers is determined by the group H of the skeleton. 
(See: Pélya, loc. cit., sections 57-59, pp. 214-17). 
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It may happen that more than one line joins point 7 to point 4, 
1.e., more than one letter is sent out in the transition j} > %. A 
simple example is given by C. Shannon and W. Weaver’s Figure 3 


(loc. cit., p. 16) in which five letters are emitted from a single 
state. Let the index 


[laste bro, ct cle = 


(where « may depend on j and &) enumerate the lines joining points 


j and &, and let es be the probability associated with each line. 
We then put 


a 
5 ad SPS (5) 
Mel 


and for the right hand sides of (5) we have again equation (4). 
The entropy of the jth state is defined as the mean value of the 
quantities —lg De for that state: 


m a 
GD) 1. CH) 
= = 1 EAP es (6) 


kul pml 


and the entropy (per symbol) of the source is then given by 


[= Vs Ens (7) 


i.e., the weighted sum of the I’s. In equation (7) P, is the proba- 
bility with which the jth state occurs. If the source is ergodic we 
use the steady-state values of the P,’s which are uniquely de- 
termined as solutions of the system 


m 
Pr= P Pip Je mth An ona ie ge tlt 
j=l 
(8) 


m 


Ne Pte shi: 


k mi 


The conditions for ‘‘ergodicity’’ are given by C. Shannon and W. 
Weaver (loc. cit., p. 17) in terms of the source-graph, and algebra- 
ically by M. Fréchet (1938, p. 112) as follows: 
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1) The equation 
A(s) = det | Dj, - 954 | = 0 (9) 


shall have no other root of modulus 1 excepts = 1.* 

2) s = 1 shall be a simple root. 

In some problems dealing with graphs it might be convenient to 
use the general formula (7). We shall, however, consider only the 
simple case in which the probabilities De do not depend on j 


(emission of one symbol not dependent on previous symbols): 


Pf = PD - (10) 


We may then associate one “‘letter’’ with each pair of values (ky 4), 
and equation (7) reduces to 


m  Q(k) 


[=~ De Ne pY Ig p , (11) 
k 


al Mel 
The matrix (P j,) now has m equal rows; its roots are 


s = 1 with multiplicity 1, 
s = 0 with multiplicity m-1; 


and the source is therefore ergodic. The source-graph may, in fact, 
be reduced to a single point with several loops ‘‘joining it to 
itself.” 

Thus equation (1) expresses the fact that in a graph G we con- 
sider a set of topologically equivalent points as a ‘‘letter’’; these 
letters or symbols are chosen independently (and with probability 
w, for the kth symbol) to form messages. 

N. Rashevsky and G. Karreman (personal communication) have 
pointed out some other possibilities of realizing a discrete source 
of information by means of a graph. One may, for example, use 
sets of topologically equivalent lines in much the same way as 
was done for the points. Rashevsky and Karreman also suggested 
introducing ‘‘compound symbols’’**, i.e., choosing simultaneously 
more than one set of topologically equivalent points as a basic 


*Adding the columns in A(s) and using (4) it is seen that unity is 
always a root of equation (9). 
**Originally called by them ‘‘higher order signals.’? The word ““sig- 


nal,’’ however, is reserved by Shannon and Weaver (Joc. cit., p. 7) to de- 
note an allowable sequence of symbols. 
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unit for the alphabet. This amounts to using an alphabet composed 
of ‘‘words’’ (combinations of letters). We must then be able to as- 
Sign a probability to each word. A simple case is the following: 
We shall assume that the probability of choosing a given word is 
determined by the product of the probabilities for the letters which 
compose it. If the pth compound symbol is obtained from the pri- 
mary symbols a, 6, c,...and only those (i.e. none of the remaining 
primary symbols s,7z, u,...is used), the probability for that ‘‘word’’ 
will be 


ioe, U,V)... (m terms) (12) 


with 


o,=1l-w, (j=1,2,...,m). (13) 


E.g. we have for three classes 1, 2, and 3 (m = 8): 


p Compound Symbol Probability P, 
1 ends O, 050, 
2 1 alone u, ¥, 5 
3 Cee.” 2 tee. 
4 Ss ees Daa Ae., 
5 1 and 2 UW, Wes 
6 1 and 3 WV, Ws, 
7 2 and 3 VU, WU, 
8 1, 2, and 3 UW, U, Us, 


In this procedure the number of different letters which can be 
chosen to give a word may be any integer from zero to m. Hence 
7 H ¢ b] 

the total number of possible compound symbols (including a ‘‘blank’’) 


will be 
y @ one (14) 


M=O 


and the entropy of the source per symbol now becomes 


gm 
a sae Pe teee os (15) 
Ppeml 


It is convenient to consider the probabilities P,, given by (12), 
in pairs, such that each pair differs only in the last factor. One of 
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these pairs would be written 
Pr = B, tt,» 


(16) 


where 8,, @ product involving (m-—1) of the quantities w,,..., 
Un» Vzysee*+9 VU» 1S the ‘same in both cases. From this law of for- 
ie . . . 
mation, using equation (13) and separating the B,’s again into 
pairs (as was done for the PEs) we derive the following equations: 

1) Since 


P?.+ Po*.= B, for allo , 


we have 
gm gm-1 
yoo) Be ee (17) 
P= o=1 


as, of course, should be the case. 
2) Furthermore: 


Po leiPee-sP2 Ve Pole pl goes (w Igw_ +, igv_). 


Hence: 
gm gm~1 
ms Po lg 2, =wesle wo le Oo D3 B, lg B, 
Puli o=l 
(18) 
er ng (uw, lg w+ %, lg v;)- 
jul 
From (15) and (18) we thus find the simple formula 
m 
Poa & = (w, lg w, + v, 1g v,). (19) 


jel 


Equation (19) could have been predicted without calculation: it 
represents the entropy of m independent ‘‘relays,’? each of them 
with two possible stable positions, of probabilities w, and 1~w j 
respectively. A. Rapoport’s (1955) formula (1) is the same as 
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equation (19) above for an assembly of neurons, when a neuron may 
be found in either one of two states, namely ‘‘on’’ and ‘‘off.”? 

In deriving equation (19) we have not used the fact that 


Dy rt, (20) 


but since (20) holds in our case, we can draw some interesting 
conclusions on the magnitude of the quantities / and 7’: 
1) Consider the function 


m 


f(#y,---,@,)= - ye x, lg @, 
fel 
_ with 0 < a, < ale 
We find the maximum value of f(a,,..., 2,) with the supple- 
mentary condition: 


m 
eae = c lg C 
for 
a ets ie “s © im Ss m 
From this we obtain: 
Tax = Em (21) 
and 
Tr ax = ™ lg m-(m—1)lg(m-—1). (22) 


For large values of m: 


uf 
Ina =lg m+ ig e+... 0(2] 


- or 


max ma 


[ioNen,] e ee V4 48 (2). (23) 


E.g. we find from (22) for m = 3, Ig m = 1.585: 
I? = 2-755 (bits per symbol) , 
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for m = 4, lg m = 2: 
1’. = 3.245, 
max 
and for m = 64, lg m= 6: 
T? 4 = $4815. 
2) We also have the inequality 
EM SO (24) 
or 
m m 
-> v,lg2,<- Dw, lg u,. (25) 
j=l jul 


The proof of (25) is based on the following properties of the 
function 


g (a) = -a@lg a+(1- 2) lg (1-2): (26) 
a) g (1-@) = -g (2) (27) 
b) g (a) > 0 for 0 <2< 1/2 (28) 
g (ey) =< Ofort/2 <a (29) 
g (0) = g (1/2) = g (1) = 0 (30) 
c) If x and y are two positive quantities such that 
e+y<1/2, (31) 
we have 
l-a~-2y>0 
ee (33) 
and 


g(z)+9(y)~g(@+y) 
=-alga+(1-—2)lg(l-2«-ylgy+(1-y)lg(i-y) 
+(a+y) lg (a+ y)-(l-«-y) lg (l-a-y) 


(33) 
(a — a”) +y(1-2¢-y) 
ge? rt 


(y-y?)+a(1-«@-2y) 


= xl 1 
(x ~ a”) 3 (y- y’) 
nar (l-2-y)+ ay 
(l-«-y) 


The third term in the last expression (33) is positive; the first 
two terms are also positive because of (32). Hence: 
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g(x) + g(y) > gle@+y). (34) 

It follows from this that if ©, Lyye+-, & are r positive quantities 
such that 

@, +, +.-.+% <1/2, (35) 


we have 


g(@,)+9(@,) +... 4+ 9(@) > 9(@,+%,+...4+2). (36) 


r 


Now put 


F (ty) Wgyee, Uy) = — ), [w,lg w,-(-u,) lg d- ul 
j=1 
(37) 


m 


=D) 9%). 


Fmt 


We may suppose that none of the w’s is equal to zero or one, since 
otherwise the corresponding terms of the sum in (37) would simply 
cancel. 

A) If all the w, *s are < 1/2 we have from (28) and (30): 


PONG sy ta 5 oe 05.) 2-0 (38) 
which is the same as (25). 


B) At most, one of the u,’s can be > 1/2. 
Assume, for instance, that 


w,>1/2. (39) 
Then 
Wo + Wat... tu, =l-w,<1/2. (40) 
Therefore, using (27) and (36) we find: 


F (45 Wao02+) Un) = -~g(1-%,)+ De g(u,) 


ju2 


j=2 


g. & d. 
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© © 
Cc D 


>© 
o 


1 , 
FIGURE 1. m= 2; wy = Wo = Vy = %Q = i> lg 2 = 1 (bits per symbol); 
[= 2, 


4 1 ADO eel 
FIGURE 2. m= 23 Wy = % = 53 =v oni le let 6 Om 0 ie 
[?.=2/7= 1,444, 


a! 3 1 1 
FIGURE 3. m= 3; wy =u, = 73 0,=0,= 55 Wy = 35 % = 55 
2 1 ‘ (nee ee 
J=-— 1 — | = 1.500: =] Sea =] = 
q le@4+ 5 le 2= 1.500; =/+ 7 le 5 + 5 le 2 = 2.6225 


In Figures 1 to 5 we give some examples for the calculation of 
I and /’, using very simple, non-directed graphs in which all points 
are assumed to be of the same kind. 

Rashevsky and Karreman further suggested treating the case in 
which a point of a given kind may exist in several different states. 

For example, if the graph is regarded as the structural formula 
of a molecule, two points occupied by different types of atoms 
(like Cl or H) may be considered of different kinds; each point or 
atom might then exist in various ‘‘states of excitation.’’ Alterna- 
tively, we could assume certain points to be of the same kind (the 
sites of atoms), but each site might exist in various states, corre- 
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A 
Q 


1 4 2 
FIGURE 4, m= 3; w, = 53 0 = 5} Wy =, = = Uy = =o 


1 4 5 4 eG 5 
/=— ] 6+ — le—= : * = a == = == 
5 le 5 85 TPs +e le Tt, le 5 = 2.664. 


1 7 5 
FIGURE 5. m=4; w,=wv,= Uv, 33 Ra SE I re Dy ire 
3 3 5 8 21 Coes 8 
oa S| Be Pe ae = Set lp as ga lp ae 
1, = 55 SE 8 5 1.549; / Rates See ee aes 2.586 


sponding to different atoms of equal valence. In Figure 2, e.g., 
the central point could be a C-atom, and each of the peripheral 
points either a Cl-atom or a H-atom (or some other monovalent 
substituent). In this way we obtain the average information content 
not of a single molecule, but of a group or class of molecules. We. 
wish to emphasize, however, that we are not, at present, trying to 
calculate the information content of any molecules, although this 
might be done subsequently. We are dealing only with graphs. 

Let each one of the n, points belonging to a topological set be 
capable of existing in 8, different states or ‘‘colors.’’ Referring to 


Figure 4, we may have, for exampie: 
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Class Points of Class Possible Colors 
s i A Red: v, = 1 
Blue: v,= 2 
Green: v, = 3 
2 Eee Purple: v, = 1 
White: v, = 9 
3 D,E Violet: v, = 1 
Yellow: v, = 2 


3 


Some of the colors allowed in different classes could be equal 
(e.g. two points of the graph may be topologically different but of 
equal degree and hence occupied by the same kind of atoms). For 
simplicity, however, we shall not consider this possibility. 

Call oP the probability that a particular point (belonging to the 
jth class) should be in the v ,th state (v, ig Pe eee '8;)- 


S: 


x weP = 1. (41) 


Viel 
There are n, topologically equivalent points in the jth class. It 
is logical to assume that the probability of choosing any one of 


them, regardless of its state, should be equal to = . Thus, ina 


random choice, the probability of selecting a particular point, in 
the v j th state, will be 


1 


De gee he: 
P= u,oo4 = 


Od (42) 


Taking only one point from the graph, the total number of pos- 
sible choices is given by n 73; for each class, and hence by 


M~)" n,'8, (43) 


a 


for the whole graph. Thus the index & in equation (42) runs from 1 
to M, and we have 


M m 


4 
Dé Prva aes : vP a1. (44) 
yas 


k =1 Vy=t 
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It is no longer possible, now, to consider the topological sets of 
points as the basic units of our alphabet, since several distinct 
arrangements may occur within one class. Rather, we shall assume 
that choosing a ‘‘letter’? means selecting a particular point in a 
given state; we have then an alphabet composed of M different 
letters. 

In this case the information content of the graph, J, becomes: 


k=l =1 View] 
: (45) 
m sj 
=Ign- de w, 3 oe yD, 
j=1 vi=1 


If each of the points can exist in one state only, equation (45) 
reduces to 


T=lgn, 


the amount of information involved in selecting one point from a 
set of n equiprobable points. Equation (45) can also be written 


[== ou, lew,~ Yo vj le > su, > Pte OP (46) 
=l 


ful j=l 1 fud 


‘ 


We may say that the first term in (46) represents the ‘‘purely 
topological’? part of 7, the second term is the entropy due to un- 
certainty of choice between the n j points of a class, and the third 
term arises from the fact that each point may be found in one of 
various possible ‘‘colors.”’ 

It is clear, however, that the idea of topological information con- 
tent loses much of its meaning once we agree to consider the ~ 
choice of a particular point as a symbol. 

We shall discuss briefly the possibility of using ‘‘compound 
symbols’? in this more general case. 

If we could arrange our M letters into words in a similar manner 
as was done before, we would have altogether 2” possible com- 
binations. Yet, two ‘‘letters’’ corresponding to the same point 
(but in different states) cannot be so combined. Choosing any 
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number r= 0, 1, 2,..., 2 Of points from the graph, each of the 
selected points being in a certain state or ‘‘color,’’ we obtain 
altogether 


Ka 7] G+) (477) 
j=l 


possible ‘‘words.’’* In fact, each point of the jth class may either 
not be chosen at all, or else chosen in any of ‘8, different states, 
j.e., altogether it may be disposed of in 8, + 1 ways. Hence equa- 
tion (47). 

We also have 


m m 
lg 2” = Ne n,8,> Ig k= )° n, lg (s,+1), 
j=i j=l 
and therefore 
Ka" (48) 


as should be the case. 


We shall, however, consider only a more restricted and very 
simple type of ‘‘word formation’’: it consists in taking simultane- 
ously all the n points of the graph, each one of them in a given 
state. The total number of words that can be formed in this way is 
reduced to 


L= J] s*i. (49) 

jet 
For the probability of a word we use again the multiplication law: 
et Pe (si=.15 8.3.35 Lys (50) 


where we have replaced the indices j (= 1, 2,..., m) and v, by p 
(now running from 1 to n) and vi (= 1, 2,...,:'8,); 6% is the same 
as oP for all points belonging to the jth class. Each of the in- 
dices vy, in (50) may vary independently from 1 to Crips 


*Or perhaps K ~ 1 words, depending on whether or not we want to in- 
clude the case of ‘*no choice’? (blank). 
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Referring again to our example above: the probability of the 
following ‘‘word’’ 


point A: red 
point B: white, point C: white 
point D: violet, point E: yellow 


It is not difficult to show that in this case we have 


= Pamal (51) 


and for the information content, now called 7’: 


n 


bel ig 


m OF 

ee (V;) (3) 
Fa, Ya aah, 
j=1 Vi= 1 


a result which is again intuitively obvious. 

In connection with equations (49) and (50) [and also (47)] the 
following should be pointed out: If within one class we had two 
arrangements which could be obtained from one another simply by 
permuting the colors on the points (e.g. point number one—red, 
point number two—green; vs. point number one—green, point num- 
ber two—red), they were regarded as different. It is equally legiti- 
mate, and perhaps even preferable, to count them as equal (‘‘one 
point red, and one point green’’). If this is done, the number of 
different arrangements, or complexions, within one class is de- 
termined only by how many points are in each state. This number, 
ry , is given by the expression 

r= Ke ye IN aged ian 
j so nj) (8,- 1)! 


L 
pee Pale P= = 
pul 


(52) 


(see, e.g., Whitworth, 1925, Proposition XXVI, page 97, or Ra- 
shevsky, 1955a, pp. 122-24). The total number of possible words 
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then becomes 


m i(n,+8,-1 
= j : 53 
nt (ue) os 
j=l j 
If there are n, , points in state 1, n, » in state 2, etc., so that 
ai} 
Ds 1; vide (j= 1, 2, , m) 
viel (54) 
Niv; 20, 


the probability of one complexion within the jth class is given by 
the product of the corresponding oP *s multiplied, however, by 
the multinomial coefficient 


n,! 


! ! i= 
ek Ad eee 
That is, each complexion has its probability represented by one 


term in the development of 
[oC eb +5... Oia 


(see, e.g., ter Haar, 1954, p. 23, or Chrystal, 1906, $12, pp. 
14-15). The probabilities of all possible words are then obtained 
by multiplication over the classes, i.e. by expanding the product 
of polynomials 


m 53 n 
TT ap 


For the information content, /’’, we find in this case: 


is: = m [oo] 77,1 [6@)]"7,2 [oSP]"i,s; 
Tet 5° |) load se aes ae 
i 
mp 


! ! ! 
shel Bs My gitisial ¢ & Bol 
(55) 
n;! 
i 
ae Peat are ! : 


* eee nN. ° 
i,2 +S; 
The inner sum in (55) refers to the jth class. It has one term 


for each set of values Mir My greeey Meg satisfying the relations 
> > J 
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(54), i.e. altogether r, terms (except that those terms for which the 
multinomial coefficient is unity vanish). 


As a simple example consider the graph of Figure 2. Assume 
that point A exists in one state only (e.g. White), but that each of 
the four points B, C, D, E may be either Red or Green, with proba- 
bility bp and 6, =1~ bp respectively. Then: 


I’ = ~—4(b, lg bp + b, lg bays 
T’ = T’ ~[4 db b¢ lg 4+6 6% 02 lg 644 5, 0% Ig 4] 
=I’ ~2 bp b, [4(bR + 82) +3 dp 6, (1+ lg 3)]. 
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ON THE DIFFUSION OF IONS IN MEMBRANES 
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A simple model of a membrane is used to obtain relations among five 
measurable quantities: the inward and outward ionic fluxes, the internal 
and external ionic concentrations, and the difference of electrical po- 
tential across the membrane. The Goldman equation is generalized to 
arbitrary geometrical shapes. 


Developments in experimental techniques within the last few 
years have made it possible to obtain reasonably accurate values 
for five parameters of basic importance for the study of ionic per- 
meability: the inward and outward ionic fluxes, the internal and 
external ionic concentrations, and the potential difference across 
the membrane. In this paper we shall derive some useful relations 
among these quantities, basing our considerations on a very simple 
model of the membrane. 

By ‘‘membrane’’ we understand a closed region of inhomogeneity 
in the ionic concentrations and electrical potential, which separates 
an ‘‘inside’’ region from an ‘‘outside’’ region. In these regions, 
denoted respectively by z and o, the concentrations and potential 
will be assumed constant. 


I. The Flux Equation 
We assume that the diffusion of ions in the membrane is described 


by the well-known equation 


\ 
Je sy (v. + a c vv) (1) 


F 
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or, as we prefer to write it, 
J=kZ"* V(cZ). (2) 


Here 


2 <r 


J = net ionic flux in mols cm™* sec 
3 


ce = ionic concentration in mols cm™ 
V =-electrostatic potential in volts 

u = ionic mobility in cm/sec per volt/cm 
2F = charge/mol of electrolyte in coulombs 


with the abbreviations 


k =-—uRT/F 


ak 
z= exp ( v). 


The boundary conditions are 


at o: C= Co, Z = Zo; 
at 2: C=-C;, Ze=1. 


(3) 


Integrating equation (2) and applying the boundary conditions, we 
find, assuming that wu is constant, 


i Byeiene eT eee) (4) 


i 


where dr indicates an element of the path of integration. We as- 
sume that a steady state obtains; in this case, integrating over the 
boundaries 


fs -dS = f J- dS = M = constant. (5) 


Put J = Mj; then equation (4) yields the flux equation 


k 
ME a ee Zo — Ci). (6) 


J Zj- dr 


u 


Il. The Flua-Ratio Equations 
The use of radioactive isotopes makes it possible to measure 
the following quantities under conditions which are equivalent ex- 


DIFFUSION OF IONS IN MEMBRANES 257 
cept for the radioactivity: 


the inward flux J;, =-J for c,;=0 
the outward flux J,,, =-J for co = 0. 


Thus in such experiments we may regard the form of Z as unchanged 
and also the form of j. Indeed, we shall see that j depended only 
on the geometry of the membrane. 
We therefore find immediately from equation (6), for a single ion, 
Min Ce 


ae (7) 


Mout Ci 
and for different ions of the same charge, 


1 1 
MSP 6 oD 


YOu > ga) pe). (8) 
in fo) 
(1) Cie CH) 

Le aes 5 

MO) ~ 42) @@)° (9) 
out t 


These relations are valid for ions moving under the influence of 
potential and concentration gradients alone and are independent of 
the geometry of the membrane. 


Ill. The Generalized Goldman Equation 

It is possible to obtain some information, quite useful in prac- 
tice (Hodgkin, 1951; Teorell, 1953; Ussing, 1953), without knowl- 
edge of the form of the potential gradient. However, in order to ob- 
tain an explicit expression for the flux, further consideration of the 
nature of the potential is evidently required. Previous treatments 
of this question by Planck (1890), Teorell (oc. cit.), and Polisser 
(1954) have involved applying our equation (1) to each ionic species 
present. Experimental observations show, however, that equation 
(1) is not sufficient to account for the movement of sodium ions 
in invertebrate nerve, although it may account for the movement of 
potassium ions. The movement of other ions has not yet been ade- 
quately investigated (see Hodgkin, Joc. cit.). 

Planck employed Poisson’s equation, where the Laplacian of V 
is proportional to the difference between the total cationic and 
anionic concentrations. Sitte (1934) argued that this difference 
would in all cases be small and that a reasonable approximation 


258 JOSEPH D. HARRIS 


might therefore be obtained by putting V V = constant. Goldman 
(1943) arrived independently at the same conclusion, basing his 
reasoning on an analogy with the copper copper-oxide rectifier. 
These investigators worked exclusively with infinite plane mem- 
branes. The flux equation derived by Goldman on this basis (loc. 
cit., see also below) has been extensively used by Hodgkin and 
his collaborators in recent years. In view of the proven usefulness 
of the Sitte-Goldman approach, it is of interest to seek a generali- 
zation of their assumption for membranes of arbitrary shape. We 
shall show that the required generalization can be obtained for the 
case V2V =0. Physically this implies an electrically neutral 
membrane with, for example, single or double layers of charge at 
the boundaries (‘‘condenser membrane’’), We regard this as a 
purely formal postulate, to be justified by the comparison of its 
implications with experiment. 

We require one further assumption, which seems reasonable phys- 
ically, namely 


curl J = 0. (10) 
A simple calculation from equation (1) yields 
curl J=-u2(VcxVV), (11) 


and hence equation (10) implies that Vc, VV and V (cZ) have the 
same direction. 
We now introduce generalized orthogonal coordinates qi, go, 93, 
(see, e.g., Margenau and Murphy, 1948) such that 
a) the surfaces g; = constant coincide with the surfaces cZ = 
constant 
b) go, gg are arbitrary functions forming a system of orthogonal 
surfaces together with q,. 


dx \? 5y\? Sz \2 
Bet.02 = [= ph ie BOE s ais 
Yh en = ee and u, = unit vector orthogo- 


nal to the surfaces g, = constant (k = 1, 2, 3). Now VV is parallel 

to V(cZ) so that VV is orthogonal to the surfaces gi = constant. 
Since 

ui OV up SV ug SV 

SS ae 

Q1 591 Q2 5q2 Qs Sq 


it follows that V = V(q;). Similarly c = c(q1). Since J = Z7! V(cZ), 
it follows in the same way that J = Ju, with J= |J|. 


(12) 
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We now have, as postulated, 


it 5 5V 
ye (220: =o (13) 
Q:1Q003 dq1 Q1 dq 
and therefore, on integration, after introducing Z 
Q1 
Z =A exp < Fi(qe, gs) dqi}. (14) 
9293 
Our steady-state assumption (5) is equivalent to 
1 5(J 
div J= See as) 5, (15) 
919293 541 
so that 
F 
Je fader ga) (16) 
Q2Qs 


Here A is a constant of integration, and F'1(qe2, 73),F'2(qe2, 93) are 
undetermined factors. 

Now since J-dr=Ju,;-dr= JQ, dqi, we are able to evaluate 
the definite integral appearing in equation (6). Thus 


Ss sf Q1 F'oQ1 
-dr=A F d d 1% 
pz8 = fot ' S003 eo. ae 


or, with equation (14), 


o F'2(q2, 93) 
.dr = —————. (Z,, - 1). 18 
if oe F'\(q2, qs) ( ee 


i 


But from equations (16) and (5) 


F ; 
Heal ty de Fo(ga> 93) Q203dqodq3, (19) 
l f Q293 


and from equation (14) and the boundary conditions (3) 


F1(qg2, 93) = aay (20) 
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Using equation (4) to eliminate F; and F'2 from equations (18), (19) 
and (20), we have 


In Z, 


M=kT On Lig eC, (21) 
ee ) 
d : 
where [* = lence a purely geometrical factor. Finally, 
fae 
29s 


expanding the abbreviations, equation (21) reads 
M r ies {= V (22) 
= as Ss eee aig (ee Ort 
Uz SF Co Xp PT 
Oxpy=— ee 
RT 


This equation shows the»same dependence on the potential differ- 
ence V, as the equation derived by Goldman. We therefore call it 
the generalized Goldman equation. 

In problems of biological interest the ratio of the membrane 
thickness d to the dimension of the cell is very small, ‘say of the 


order of 107*. In this case, since d -{ Qidq, we may put 
i 


ep re (23) 
QoQ QoQ rail 
approximately. Then 
Dee f Q22s Biter (24) 
d l 71249100) ad 
where A is the area of the membrane. Now 
os =-Jay = z J-ds (25) 
A A [ 


so that to a good approximation 


wz Vo oF 
Say =— es a FE ee Co exp RT Vo — C; (26) 
ye sat 3 
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A 


Goldman’s equation differs from equation (26) only in that he did 
not find it necessary to average J over the membrane surface, 
since he considered only infinite plane membranes. 


It is a pleasure to thank Professor L. J. Mullins for many stim- 
ulating and enlightening discussions. The author is also indebted 
to Professor D. ter Haar for a most valuable suggestion. 
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ERRATA 


1. To the paper by D. G. O'Sullivan entitled ‘‘Diffusion and 
Simultaneous Chemical Reactions: I. A Method for Solving the 
Equations of Some Systems in which a Fixed Concentration Exists 
at a Boundary”’ (Vol. 17, pp. 141-58). 

In equation (5) on page 144 the second equality sign should be 
replaced by a minus sign, and in the penultimate line on page 151, 
a8 should be replaced by af in the denominator of the expression 
for y. 


2. To the paper by D. G. O'Sullivan entitled ‘‘Diffusion and 
Simultaneous Chemical Reactions: I]. The Equations of Those 
Systems in which Transport Occurs from One Region to an Adjoin- 
ing Region’’ (Vol. 17, pp. 243-55). 

In the second equation that appears in Section IV of this paper 
on page 246,—A, should be replaced by + A, , and the lower limits 
of integration that have not been printed in the last equation on 
page 246 and in the penultimate equation on page 247 should both 
be zero. The final displayed mathematical expression on page 254 
should be disregarded. This is rediscussed in Part II of this series 
entitled ‘‘Diffusion and Simultaneous Chemical Reactions: II. The 
Degree of Localization Achieved in Cytochemical Staining Proce- 
dures’’ (Vol. 18, pp. 200-201). 


